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Abstract 


This  paper  reviews  some  recent:  developments  in  the  convex  analysis 
approach  to  integer  programmi ng . These  developments  are  based  on  viewing 
integer  programs  as  disjunctive  programs,  i.e.,  linear  programs  with 
disjunctive  constraints,  an  approach  which  seems  to  be  fruitful  for  0-1 
prog, ramming  both  theoretic;!  1 1 y and  practically.  On  the  theoretical  side, 
it  provides  structural  characterizations  which  offer  new  insights.  On  the 
practical  side,  it  produces  a variety  of  cutting  planes  with  desirable 
properties  and  offers  new  ways  of  combining  cutting  planes  with  enumerative 
techniques . 
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DI5 JUNCTIVE  PROGRAMMING 
by 

Egon  Balas 

,i 


1.  Introduction 

This  paper  reviews  some  recent  developments  in  the  convex  analysis 
approach  to  integer  programming.  A product  of  the  last  four  years,  these 
developments  are  based  on  viewing  integer  programs  as  disjunctive  programs, 
i.e.,  linear  programs  with  disjunctive  constraints.  Apart  from  the  fact 
that  this  is  the  most  natural  and  straightforward  way  of  stating  many 
problems  involving  logical  conditions  (dichotomies , implications,  etc.), 
the  disjunctive  programming  approach  seems  to  be  fruitful  for  zero-one 
programming  both  theoretically  and  practically.  On  the  theoretical  side, 
it  provides  neat  structural  characterizations  which  offer  new  insights. 

On  the  practical  side,  it  produces  a variety  of  cutting  planes  with  desirable 
properties,  and  offers  several  ways  of  combining  cutting  planes  with  branch 
and  bound. 

The  line  of  research  that  has  led  to  the  disjunctive  programming  approach 
originates  in  the  work  on  intersection  or  convexity  cuts  by  Young  [40],  Balas  [2], 
Glover  [23],  Owen  [36]  and  others  (see  also  [13],  [24],  [42]).  This  geometrically 
motivated  work  can  be  described  in  terms  of  intersecting  the  edges  of  the  cone 
originating  at  the  linear  programming  optimum  x with  the  boundary  of  some  convex 
set  S,  whose  interior  contains  x but  no  feasible  integer  points,  and  using  the 
hyperplane  defined  by  the  intersection  points  as  a cutting  plane.  An  early 
forerunner  of  this  kind  of  approach  was  a paper  by  Hoang  Tuy  [29], 


In  the  early  70's,  research  on  intersection  or  convexity  cuts  was 
pursued  in  two  main  directions.  One,  typified  by  [23],  [18],  [4],  was 
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aimed  at  obtaining  stronger  cuts  by  including  into  S some  explicitly  or 
implicitly  enumerated  feasible  integer  points.  The  other,  initiated  by  [3], 
brought  into  play  polarity,  in  particular  outer  polars  (i.e.,  polars  of  the 
feasible  set,  scaled  so  as  to  contain  all  feasible  0-1  points  in  their 
boundary),  and  related  concepts  of  convex  analysis,  like  maximal  convex 
extensions,  support  and  gauge  functions,  etc.  (see  also  [5  ],  [15], 

[19]).  Besides  cutting  planes,  this  second  direction  has  also  produced 
(see  [5  ],  [6  ])  a "constraint  activating"  method  (computationally  untested 
to  date)  based  on  the  idea  of  "burying"  the  feasible  set  into  the  outer  polar 
(without  using  cuts),  by  activating  the  problem  constraints  one  by  one,  as 
needed.  This  research  has  yielded  certain  insights  and  produced  reasonably 
strong  cutting  planes;  but  those  procedures  that  were  implemented  (and,  by  the 
way,  very  few  were)  turned  out  to  be  computationally  too  expensive  to  be 
practically  useful. 

In  1973  Glover  [25]  discovered  that  intersection  cuts  derived  from  a 
convex  polyhedron  S can  be  strengthened  by  rotating  the  facets  of  S in 
certain  ways,  a procedure  he  called  polyhedral  annexation.  This  was  an 
important  step  toward  the  development  of  the  techniques  discussed  in  this 
paper.  The  same  conclusions  were  reached  independently  (and  concomitantly) 
in  a somewhat  different  context  by  Balas  [7  ].  The  new  context  was  given 
by  the  recognition  that  intersection  cuts  could  be  viewed  as  derived  from  a 
disjunction.  Indeed,  requiring  that  no  feasible  integer  point  be  contained 
in  the  interior  of  S,  is  the  same  as  requiring  every  feasible  integer  point 
to  satisfy  at  least  one  of  the  inequalities  whose  complements  define  the 
supporting  halfspaces  of  S.  This  seemingly  innocuous  change  of  perspective 
proved  to  be  extremely  fruitful.  For  me  thing  it  led  naturally  and 
immediately  to  the  consideration  of  disjunctive  programs  in  their  generality 
[ 8 ],  [ 9 ],  and  hence  to  a characterization  of  all  the  valid  inequalities 


for  an  Integer  program.  By  the  same  token,  it  offered  the  new  possibility 
of  generating  cuts  specially  tailored  for  problems  with  a given  structure. 
Besides,  it  offered  a unified  theoretical  perspective  on  cutting  planes 
and  enumeration,  as  well  as  practical  ways  of  combining  the  two  approaches. 
Finally,  it  vastly  simplified  the  proofs  of  many  earlier  results  and  opened 
the  way  to  the  subsequent  developments  to  be  discussed  below. 

Besides  the  above  antecedents  of  the  disjunctive  programming  approach 
there  have  been  a few  other  papers  concerned  with  linear  (or  nonlinear) 
programs  with  disjunctive  constraints  [26],  [27].  The  paper  by  Owen  [36] 
deserves  special  mention,  as  the  first  occurrence  of  a cut  with  coefficients  of 
different  signs.  However,  these  efforts  were  focused  on  special  cases. 

The  main  conceptual  tool  used  in  studying  the  structural  properties 

of  disjunctive  programs  is  polarity.  Besides  the  classical  polar  sets, 

we  use  a natural  generalization  of  the  latter  which  we  call  reverse  polar 

sets  [10],  This  connects  our  research  to  the  work  of  Fulkerson  [20],  [21], 

[22],  whose  blocking  and  antiblocking  polyhedra  are  close  relatives  of  the 

polar  and  reverse  polar  sets.  There  are  also  connections  with  the  work  of 

Tind  [39]  and  Araoz  [ 1 ].  One  crucial  difference  in  the  way  polars  and 

reverse  polars  of  a polyhedron  are  used  in  our  work,  versus  that  of  the  above 

y/ 

mentioned  authors,  is  the  fact  that  we  "dualize"  the  reverse  polar  S'  of  a 
(disjunctive)  set  S,  by  representing  it  in  terms  of  the  inequalities  (of  the 
disjunction)  defining  S,  rather  than  in  terms  of  the  points  of  S.  It  is 
precisely  this  element  which  leads  to  a linear  programming  characterization 
of  the  convex  hull  of  feasible  points  of  a disjunctive  program. 

Except  for  the  specific  applications  described  in  sections  7 and  8, 
which  are  new,  the  results  reviewed  here  are  from  [ 8 ],  [ 9 ],  [10],  [14]. 
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We  tried  to  make  this  review  self-contained,  giving  complete  proofs  for 
most  of  the  statements.  Most  of  the  results  are  illustrated  on  numerical 
examples.  For  further  details  and  related  developments  the  reader  is 
referred  to  the  above  papers,  as  well  as  those  of  Glover  [25], 

Jeroslow  [32],  and  Balas  [11].  Related  theoretical  developments  are  also 
to  be  found  in  Jeroslow  [30],  [31],  Blair  and  Jeroslow  [17],Zemel  [41], 
while  some  related  computational  procedures  are  discussed  in  Balas  [12], 

This  paper  is  organized  as  follows.  Section  2 introduces  some  basic 
concepts  and  terminology,  and  discusses  ways  of  formulating  disjunctive 
programs  (DP).  Section  3 contains  the  basic  characterization  of  the  family 
of  valid  inequalities  for  a DP,  and  discusses  some  of  the  implications  of 
this  general  principle  for  deriving  cutting  planes.  Section  4 extends  the 
duality  theorem  of  linear  programming  to  DP.  Section  5 discusses  the  baste 
properties  of  reverse  polars  and  uses  them  to  characterize  the  convex  hull 
of  a disjunctive  set  (the  feasible  set  of  a DP).  It  shows  how  the  facets  of 
the  convex  hull  of  the  feasible  set  of  a DP  in  n variables,  defined  by  a dis- 
junction with  q terms,  can  be  obtained  from  a linear  program  with  q(n  + 1) 
constraints . 

In  section  6 we  address  the  intriguing  question  of  whether  the  convex 
hull  of  a disjunctive  set  can  be  generated  sequentially,  by  imposing  one 
by  one  the  disjunctions  of  the  conjunctive  normal  form,  and  producing  at 
each  step  the  convex  hull  of  a disjunctive  set  with  one  elementary  disjunction. 
We  answer  the  question  in  the  negative  for  the  case  of  a general  DP  or  a 
general  integer  program,  but  in  the  positive  for  a class  of  DP  called  facial, 
which  subsumes  the  general  (pure  or  mixed)  zero-one  program,  the  linear 
complementarity  problem,  the  nonconvex  (linearly  constrained)  quadratic 


program,  etc 
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The  first  6 sections  treat  the  integrality  constraints  of  an  (otherwise) 
linear  program  as  disjunctions.  When  it  comes  to  generating  cutting  planes 
from  a particular  noninteger  simplex  tableau,  the  disjunctions  that  can  be 
used  effectively  are  those  involving  the  basic  variables.  Section  7 discusses 
a principle  for  strengthening  cutting  planes  derived  from  a disjunction,  by 
using  the  integrality  conditions  (if  any)  on  the  nonbasic  variables.  Section  8 
illustrates  on  the  case  of  multiple  choice  constraints  how  the  procedures  of 
sections  3 and  7 can  take  advantage  of  problem  structure.  Section  9 discusses 
some  ways  in  which  disjunctive  programming  can  be  used  to  combine  branch  and 
bound  with  cutting  planes.  In  particular,  if  LP^,  keQ,  are  the  subproblems 
associated  with  the  active  nodes  of  the  search  tree  at  a given  stage  of  a 
branch  and  bound  process  applied  to  a mixed  integer  program  P,  it  is  shown 
that  a cut  can  be  derived  from  the  cost  rows  of  the  simplex  tableaux  associated 
with  the  problems  LP^,  which  provides  a bound  on  the  value  of  P,  often 
consideiably  better  than  the  one  available  from  the  branch  and  bound  process. 
Finally,  section  10  deals  with  techniques  for  deriving  disjunctions  from 
conditional  bounds.  Cutting  planes  obtained  from  such  disjunctions  have  been 
used  with  good  results  on  large  sparse  set  covering  problems. 




-6- 

Llnear  Programs  with  Logical  Conditions 

By  disjunctive  programming  we  mean  linear  programming  with  disjunctive 
constraints.  Integer  programs  (pure  and  mixed)  and  a host  of  other  non- 
convex  programming  problems  (the  linear  complementarity  problem,  the  general 
quadratic  program,  separable  nonlinear  programs,  bimatrix  games,  etc.)  can 
be  stated  as  linear  programs  with  logical  conditions.  In  the  present  context 
"logical  conditions"  means  statements  about  linear  inequalities,  involving 
the  operations  "and"  (A,  conjunction--sometimes  denoted  by  juxtaposition), 

"or"  (V,  disjunction),  "complement  of"  (-1,  negation).  The  operation 
"if. ..then"  (=»  , implication)  is  known  to  be  equivalent  to  a disjunction. 

The  negation  of  a conjunctive  set  of  inequalities  is  a disjunction  whose  terms 
are  the  same  inequalities.  The  operation  of  conjunction  applied  to  linear 
inequalities  gives  rise  to  (convex)  polyhedral  sets.  The  disjunctions  are 
thus  the  crucial  elements  of  a logical  condition  (the  ones  that  make  the 
const  aint  set  nonconvex) , and  that  is  why  we  call  this  type  of  problem  a 
disjunctive  program. 

A disjunctive  program  (DP)  is  then  a problem  of  the  form 
min  {cx  j Ax  > ag,  x > 0,  xel] 

where  A is  an  m x n matrix,  aQ  is  an  m-vector,  and  L is  a set  of  logical 
conditions.  Since  the  latter  can  be  expressed  in  many  ways,  there  are 
many  different  forms  in  which  a DP  can  be  stated.  Two  of  them  are 
fundamental. 

The  constraint  set  of  a DP  (and  the  DP  itself)  is  said  to  be  in 
disjunctive  normal  form  if  it  is  defined  by  a disjunction  whose  terms  do 
not  contain  further  disjunctions;  and  in  conjunctive  normal  form,  if  it  is 
defined  by  a conjunction  whose  terms  do  not  contain  further  conjunctions. 


The  disjunctive  normal  form  is  then 

(2.1) 

while  the  conjunctive  normal  form  is 


V 

heQ 


,h  ^ h 1 
A x > 

x > 0 


(2.2) 


or,  alternatively. 


>a0 

x > 0 

V (dSc^d  ), 

ieQj 


j e S 


1 Ax  > 

n 

- 

(2.2 ')  | 

f — Q 

i ° i 

1 A 

V (dXx>d  ), 

Lia>i  j 

A ...  A 

v (d1*  > d > 

Ksl  J 

Here  each  d is  an  n-vector  and  each  d^  a scalar,  while  the  sets  Q and 
Qj*  j eS,  may  or  may  not  be  finite.  The  connection  between  the  two  forms  is 
that  each  term  of  the  disjunction  (2.1)  has,  besides  the  m + n inequalities 
of  the  system  Ax  > a^,  x > 0,  precisely  one  inequality  d*x  > d^,  i eQ^, 
from  each  disjunction  j e S of  (2.2),  and  that  all  distinct  systems 
A^x  > aQ>  x > 0 with  this  property  are  present  among  the  terms  of  (2.1);  so 
that,  if  Q (and  hence  each  Q.,  j eS)  is  finite,  then  |q  | = | tt  Q.|»  where 

J jes  J 

rr  stands  for  cartesian  product.  Since  the  operations  A and  V are  distributive 
with  respect  to  each  other  [i.e.,  if  A,  B,  C are  inequalities,  A A (B  V C)  = 

= AB  V AC,  and  A v (3C)  = (A  V B)  A (A  V C)],  any  logical  condition  involving 
these  operations  can  be  brought  to  any  of  the  two  fundamental  forms,  and  each 
of  the  latter  can  be  obtained  from  the  other  one. 

We  illustrate  the  meaning  of  these  two  forms  on  the  case  when  the  DP  is 
a zero-one  program  in  n variables.  Then  the  disjunctive  normal  form  (2.1)  is 

V (Ax  > 0,  x > 0,  x «*  xh) 
heQ 


J 


when.  x\...,x*^'  is  the  set  of  all  0-1  points,  and  |q|  » 2n;  whereas  the 
conjunctive  normal  form  is 

Ax  > 0,  x > 0,  (Xj  = 0)  V (Xj  = 1),  j = l,...,n. 

Once  the  inequalities  occurring  in  the  conjunctions  and/or  disjunctions 
of  a DP  are  given,  the  disjunctive  and  conjunctive  normal  forms  are  unique. 
It  is  a fact  of  crucial  practical  importance,  however,  that  the  inequalities 
expressing  the  conditions  of  a given  problem  can  be  chosen  in  more  than  one 
way.  For  instance,  the  constraint  set 


3x. 

+ 

x„  - 

2x.  + x. 

< 

1 

2 

3 4 

x. 

+ 

x.  + 

x_  + X, 

< 

1 

2 

3 4 

XJ 

= 

0 or 

1,  j - 1, 

• • • > 

whs:.-!  put  in  disjunctive  normal  form,  becomes  a disjunction  with  2 = 16  terms 

but  the  same  constraint  set  can  also  be  expressed  as 


3x,  + x0  - 2x„  + x.  < 1 

1 2 3 4 — 


4 

V 

i=l 


Xj  - 0,  j + i j 


0,  Vj) 


which  gives  rise  to  a disjunction  with  only  5 terms 
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A constraint  B is  said  to  be  a consequence  of,  or  Implied  by,  a constraint 

A,  if  every  x that  satisfies  A also  satisfies  B.  We  are  interested  in  the 

family  of  inequalities  implied  by  the  constraint  set  of  a general  disjunctive 

program  (DP).  All  valid  cutting  planes  for  a DP  belong  of  course  to  this 

family.  On  the  other  hand,  the  set  of  points  satisfying  all  the  inequalities 

in  the  family  is  precisely  the  convex  hull  of  the  set  of  feasible  solutions  to 

the  DP  • A characterization  of  this  family  is  given  in  the  next  theorem, 

which  is  an  easy  but  important  generalization  of  a classical  result.  Let 
n n h mhX  n h ’"h 

x e R , arSR  , ar^  e R,  A eR  , a^  eR  , h cQ  (not  necessarily  finite)  and 

let  a^  be  the  j-th  column  of  A , h eQ,  jcN  = {1, . . . ,n}. 

Theorem  3.1.  The  inequality  a x > ^ is  a consequence  of  the  constraint 


h ti  ^ 

if  and  only  if  there  exists  a set  of  6 eR  ,0  >0,  heQ,  satisfying 


(3.2)  a > 6hAh  and  aQ  < 0^  , VheQ*. 

* h 

where  Q is  the  set  of  those  heQ  such  that  the  system  A x > a^,  x > 0 is 

consistent. 


Proof,  a x > Oq  is  a consequence  of  (3.1)  if  and  only  if  it  is  a 

* 

consequence  of  each  term  h€Q  of  (3.1).  But  according  to  a classical  result 
on  linear  inequalities  (see,  for  instance.  Theorem  1.4.4  of  [37],  or 
Theorem  22.3  of  [36]),  this  is  the  case  if  and  only  if  the  conditions  stated 
in  the  theorem  hold.  || 

•n* 

Remark  3.1.1.  If  the  i-th  inequality  of  a system  heQ  of  (3.1)  is 
replaced  by  an  equation,  the  i-th  component  of  0 is  to  be  made  unconstrained. 
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If  the  variable  Xj  In  (3.1)  is  allowed  to  be  unconstrained,  the  j-th 

h h ^ 

inequality  of  each  system  a > 0 A , h*Q  , is  to  be  replaced  by  the  corre- 
sponding equation  in  the  "if"  part  of  the  statement. 

With  these  changes.  Theorem  3.1  remains  true. 

An  alternative  way  of  stating  (3.2)  is 


(3.3) 


> sup*  ehfll! » 

J hCQ  J 

“o  - ehflo  * 

hcQ 


J « N 


Since  Q c Q,  the  if  part  of  the  Theorem  remains  of  course 
* 

valid  if  Q is  replaced  by  Q. 

Since  (3.3)  defines  all  the  valid  inequalities  for  (3.1),  every 
valid  cutting  plane  for  a disjunctive  program  can  be  obtained  from  (3.3) 
by  choosing  suitable  multipliers  0^  . If  we  think  of  (3.1)  as  being  expressed 
in  terms  of  the  nonbasic  variables  of  a basic  optimal  solution  to  the  linear 
program  associated  with  a DP,  then  a valid  inequality  qx  > c*q  cuts  off  the 
optimal  linear  programming  solution  (corresponding  to  x^  = 0,  JcN)  if  and 
only  if  qtq  > 0;  hence  c*q  will  have  to  be  fixed  at  a positive  value.  Inequalities 
with  cxq  < 0 may  still  cut  off  parts  of  the  linear  programming  feasible  set,  but 
not  the  optimal  solution  x = 0. 

h h 

The  special  case  when  each  system  A x > a^  , h eQ,  consists  of  the 
single  inequality  a^x  > a^  (a*1  a vector,  a^  a positive  scalar)  deserves 
special  mention.  In  this  case,  choosing  multipliers  0^  = 1/a^Q,  heQ,  we 
obtain  the  inequality 


(3.4) 


E (max  a'j/a.  ) x > 1 , 

jeJ  heQ  J U J 
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r 


which  (for  Q finite)  is  Owen's  cut  [36]  . it  can  also  be  viewed  as  a 
slightly  Improved  version  of  the  intersection  cut  from  the  convex  set 

S - {«  | «hx  < «h0,  h .Q)  , 

which  has  the  same  coefficients  as  (3.4)  except  for  those  (if  any)  j e J such 
that  < 0,  Vh  cQ.  For  the  latter,  the  intersection  cut  from  S has  zero 
coefficients  whereas  the  corresponding  coefficients  of  (3.4)  are  negative. 

Whenever  all  the  coefficients  of  (3.4)  are  positive  (in  terms  of 

Intersection  cuts,  this  corresponds  to  the  case  when  S is  bounded),  (3.4) 

is  the  strongest  inequality  implied  by  the  disjunction  V (a^x  > a,  _);  in 

hCQ  “ “° 

the  presence  of  negative  coefficients,  however,  (3.4)  can  sometimes  be 
further  strengthened. 

Due  to  the  generality  of  the  family  of  inequalities  defined  by  (3.3), 
not  only  can  the  earlier  cuts  of  the  literature  be  easily  recovered  by  an 
appropriate  choice  of  the  multipliers  8*1  (see  [8]  for  details), 
but  putting  them  in  the  form  (3.3)  indicates,  by  the  same  token,  ways  in 
which  they  can  be  strengthened  by  appropriate  changes  in  the  multipliers. 

A significant  new  feature  of  the  cutting  planes  defined  by  (3.3) 
consists  in  the  fact  that  they  can  have  coefficients  of  different  signs. 

The  classical  cutting  planes,  as  well  as  the  early  intersection/convexity 
cuts  and  the  group  theoretic  cutting  planes  (including  those  corresponding  to 
facets  of  the  corner  polyhedron),  are  all  restricted  to  positive  coefficients 
(when  stated  in  the  form  > , in  term?  of  the  nonbasic  variables  of  the  tableau 
from  which  they  were  derived).  This  important  limitation,  which  tends  to 
produce  degeneracy  in  dual  cutting  plam\  algorithms,  can  often  be  overcome 
in  the  case  of  the  cuts  obtained  from  (3.3)  by  an  appropriate  choice  of  multi- 
pliers.  ( 
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Anoth«r  important  feature  of  the  principle  expressed  in  Theorem  3.1 
for  generating  cutting  planes  is  the  fact  that  in  formulating  a given  Integer 
program  as  a disjunctive  program,  one  can  take  advantage  of  any  particular 
structure  the  problem  may  have.  In  section  7 we  will  Illustrate  this  on 
some  frequently  occurring  structures.  We  finish  this  section  by  an  example 
of  a cut  for  a general  mixed  Integer  program. 


Example  3.1.  Consider  the  mixed  integer  program  whose  constraint 
set  is 


Xj  - .2  + ,4(-x3)  + 1.3(-x4)  - .01(-x5)  + .07(-x6) 
x2  - .9  - .3(-x3)  + .4(-x4)  - .04(-x5)  + .1  (-x6) 

Xj  > 0,  J - 1, . . . ,6,  Xj  Integer,  j » 1 4. 

This  problem  is  taken  from  Johnson's  paper  [35],  which  also  lists  six 
cutting  planes  derived  from  the  extreme  valid  inequalities  for  the  associated 
group  problem: 


75  x3  + .875x4  + 
778x3  + .444x^  + 
333x,  + ,667x4  + 
50  x3  + x4  + 
444x3  + .333x4  + 
394x3  + ,636x4  + 


.0125x3  + .35  Xg 
.40  x3  + .lllx6 
.033x3  + .35  Xg 
.40  x^  + .25  Xg 
,055x3  + ,478xg 
.346x3  + . 155xg 


> 1 
> 1 
> 1 
> 1 
> 1 
> 1 . 


The  first  two  of  these  inequalities  are  the  mixed-integer  Gomory  cuts 
derived  from  the  row  of  and  x^  respectively.  To  show  how  they  can  be 
improved,  we  first  derive  them  as  they  are.  TO  do  this,  for  a row  of  the 
form 


x 


1 


E 

Jej 


aU( 


■Xj) 


n 


. 'M(  » MS  U;  & 


•*  S L 
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with  integer-constrained  for  JcJ^,  continuous  for  J c Jj,  one  defines 
*ij  " *lj"^*lj^*  ^ * JU{0],  end  <PiQ  - f^^, 


£ij 

J«jJ  - £J « Jx  I 

1 £10  - £ljJ 

fiJ  ■ 1 

J * - (J  « Jj  | 

1 £i0  < fijl 

•ij 

J • Jj  * 

Then  every  x which  satisfies  the  above  equation  and  the  integrality 
constraints  on  Xj,  jeJ^Ufi),  also  satisfies  the  condition 

yt  “ <P10  + ^E  ^jC-Xj)  , yt  integer  . 

For  the  two  equations  of  the  example,  the  resulting  conditions  are 

yx  - .2  - ^(-Xj)  - .7(-x4)  - .Ol(-Xg)  + .07(-x6),  y^  integer 

y2  - .9  + ^(-Xj)  + .4(-x4)  - ,04(-x5)  + .1  (-xg),  y2  Integer. 

Since  each  y^  Is  Integer- cons trained,  they  have  to  satisfy  the 
disjunction  y^  < 0 V > 1.  Applying  the  above  theorem  with  multipliers 
8^  • 1/s^Q  then  gives  the  cut  (3.4)  which  in  the  two  cases  1 ■ 1,2  is 

.6  , . 7 . . 01  , . 07  v . 

.8  *3  .8  x4  + .8  x5  + .2  x6  ^ 1 

and 

«7  . .4  i .04  , »1  ^ i 

.9  3 .9  X4  .1  5 .9  X6  - 1 * 

These  are  precisely  the  first  two  inequalities  of  the  above  list. 

Since  all  cuts  discussed  here  are  stated  in  the  form  > 1,  the  smaller 

the  j-th  coefficient,  the  stronger  is  the  cut  in  the  direction  j.  We  would 

thus  like  to  reduce  the  size  of  the  coefficients  as  much  as  possible. 

Now  suppose  that  instead  of  y^  < 0 V y^  > 1,  we  use  the  disjunction 

1,1 5 0)  v fc ! o} 


which  of  course  is  also  satisfied  by  every  feasible  x. 


It.  .«  ft 


- i.'t- 


T^en,  applying  Theorem  3.1  with  multipliers  5,  5 and  15  for  < 0, 
y^  > 1 and  > 0 respectively,  we  obtain  the  cut  whose  coefficients 
are 

raax  S 5x  (-._6i  5,x.,6+15x.Cr.41  I . .3 

“ 1 5X.2  ’ 5 X .8+  15  x (-.2)  J 

f 5 x (-.7)  5 x .7  + 15  x (-1.3)  ) 

"ax  i 5 x .2  ' 5x.8+l5x(-.2)  ; "3-5 

i 5 X-.2  • 5x.8+X5x(-.2)  ; ‘2 

SP.f-SH.  Sx(-.07>f  UiM?)  j . 35 

max  l 5 x »2  ’ 5 x .8+15  x (-.2)  J * ’ 


that  is 


-3x„  - 3.5x.  + .2x,  + .35x,  > 1 . 

3 4 5 6 — 


The  sum  of  coefficients  on  the  left  hand  side  has  been  reduced  from 
1.9875  to  -5.95. 

Similarly,  for  the  second  cut,  if  instead  of  < 0 V y^  > 1 we  use 


the  disjunction 


y2  < o 

xx  > 0 


v (y,  > i}  » 


with  multipliers  10,  40  and  10  for  y2  < °»  xi  > 0 and  y2  > 1 respectively. 


we  obtain  the  cut 


-7X3  - 4x^  + .4x<-  - Xg  > 1 . 

Here  the  sum  of  left  hand  side  coefficients  has  been  reduced  from 
1.733  to  -11.6. 
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4.  Duality 

In  this  section  we  state  a duality  theorem  for  disjunctive  programs, 
which  generalizes  to  this  class  of  problems  the  duality  theorem  of  linear 
programming. 

Consider  the  disjunctive  program 

I - min  ex 
o 


| Ahx  > bh 


<P) 


V 

hcQ  | * > 0 


where  Ah  le  • matrix  and  bh  a vector,  ¥ hcQ. 

We  define  the  dual  of  (P)  to  be  the  problem 


w - max  w 
o 


<D) 


A 

hcQ 


w - uhbh  < 0 

uhAh  < c 

uh  > 0 


Th.  cn.tr.lnt  .«  of  <»>  require.  ««•*  »\  heq.  to  s.tl.fy  the  corre.pondlng 
bracketed  system,  and  w to  satisfy  each  of  them. 
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We  will  assume  the  following 
legularltv  condition: 


(Q*  i 0,  Q\Q**  * 0)  • Q*\Q**  4 0; 

l.e..  If  (P)  Is  feasible  and  (D)  is  infeasible,  then  there  exists  heQ  such  that 

\ i 0,  uh  - 0. 

Theorem  4.1.  Assume  that  (P)  and  (D)  satisfy  the  regularity  condition. 
Then  exactly  one  of  the  following  two  situations  holds. 

1.  Both  problems  are  feasible;  each  has  an  optimal  solution  and 

z ■ w . 
o o 

2.  One  of  the  problems  Is  infeasible;  the  other  one  either  is 
Infeasible  or  has  no  finite  optimum. 

Proof.  (1)  Assume  that  both  (P)  and  (D)  are  feasible.  If  (P)  has  no 

finite  minimum,  then  there  exists  hcQ  such  that  X 4 0 and  x e X^  such 

h " 

that  cx  < 0.  But  then  ■ 0,  l.e.,  (D)  is  infeasible;  a contradiction. 

Thus  (P)  has  an  optimal  solution,  say  x.  Then  the  Inequality  cx  > z 

o 

Is  a consequence  of  the  constraint  set  of  (P) ; l.e.,  x < X^  implies  cx  > zq, 

¥ heQ.  But  then  for  all  heQ*,  there  exists  u^eU.  such  that  u^b^  > z . 

h — o 

Further,  since  (D)  Is  feasible,  for  each  heQ\Q*  there  exists  u^eU^;  and 
since  = 0 (for  h«Q\Q*),  there  also  exists  u*1^  such  that  u**bh  > 0, 

¥ heQ\Q*.  But  then,  defining 

uh(X)  - u1  + X^*,  h«Q\Q*  , 

for  X sufficiently  large,  u^(X)«Uh,  u^(X)b^  > zq,  ¥ heQ\Q*. 

Hence  for  all  heQ,  there  exist  vectors  u*1  satisfying  the  constraints 
of  (D)  for  w ■ zq.  To  show  that  this  is  the  maximal  value  of  w,  we  note 
chat  since  x 1"  optimal  for  (P),  there  exists  heQ  such  that 


L 
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cx  «=  min{cx|x  e X^}  . 

But  then  by  linear  programming  duality, 
ex  * n»ax(u''b''|u''eUji) 

■ max{w|w  - u^b*1  < 0,  u^eU^} 

> max{w|  A (w  - u^b*1  < 0,  u^cU.  )} 
hcQ 

l.e.,  w < z , and  hence  the  maximum  value  of  w Is  w = z . 

— o o o 

(11)  Assume  that  at  least  one  of  (P)  and  (D)  is  Infeasible.  If  (P) 

Is  Infeasible,  X^  = 0,  ¥ hcQ;  hence  for  all  hcQ,  there  exists  u^eU^  such 

that  u^bh  > 0. 

If  (D)  Is  Infeasible,  we  are  done.  Otherwise,  for  each  hcQ  there 
exists  ueU^.  But  then  defining 

uh(A)  = uh  + Xu''  , heQ, 


h h h 

u (X.)eU^,  heQ,  for  all  X ">  0,  and  since  u b > 0,  ¥ heQ,  w can  be  made 
arbitrarily  large  by  Increasing  X;  l.e.,  (D)  has  no  finite  optimum. 

Conversely,  If  (D)  Is  Infeasible,  then  either  (P)  Is  Infeasible  and 
we  are  done,  or  else,  from  the  regularity  condition,  Q*\Q**  ^ 0;  and  for 
heQ*\Q**  there  exists  x e X^  and  x e X^  suih  that  cx  < 0.  But  then 


x(p)  = x + px 

Is  a feasible  solution  to  (P)  for  any  p > 0,  and  since  cx  < 0,  z can  be 
made  arbitrarily  small  by  increasing  p;  l.e.,  (P)  has  no  finite  optimum. 

Q.E.D. 
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* 

The  above  theorem  asserts  that  either  situation  l or  situation  ? hold* 

• . for  (P)  and  (D)  lf^  the  regularity  condition  Is  satisfied.  The  following 

Corollary  shows  that  the  condition  Is  not  only  sufficient  but  also  necessary. 

Corollary  4,1,1.  If  the  regularity  condition  does  not  hold,  then 
If  (P)  Is  feasible  and  (D)  is  infeasible,  (P)  has  a finite  minimum  (i.e., 
there  is  a "duality  gap"). 

Proof.  Let  (P)  be  feasible,  (D)  infeasible,  and  q*\q**  * 0,  l.e., 
for  every  heQ*,  let  / 0.  Then  for  each  heO*,  min{cx|x  c X^}  Is  finite, 
hence  (P)  has  a finite  minimum.  Q.E.D. 

Remark.  The  theorem  remains  true  If  some  of  the  variables  of  (P) 

[of  (D) ] are  unconstrained,  and  the  corresponding  constraints  of  (D)  [of  (P)] 
are  equalities. 

The  regularity  condition  can  be  expected  to  hold  In  all  but  some 
rather  peculiar  situations.  In  linear  programming  duality,  the  case  when 
both  the  primal  and  the  dual  problem  Is  Infeasible  only  occurs  for  problems 
whose  coefficient  matrix  A has  the  rather  special  property  that  there  exists 

I 

x { 0,  u t 0,  satisfying  the  homogeneous  system 
Ax  > 0 , x > 0 
uA  < 0 , u > 0 

In  this  context,  our  regularity  condition  requires  that,  if  the  primal  problem 
Is  feasible  and  the  dual  Is  infeasible,  then  at  least  one  of  the  matrices 

A^*  whose  associated  sets  D,  are  infeasible,  should  not  have  the  above 

n 

mentioned  special  property. 

Though  most  problems  satisfy  this  requirement,  nevertheless  there 
are  situations  when  the  regularity  condition  breaks  down,  as  Illustrated  by 


the  following  example. 


Consider  the  disjunctive  program 


min  -Xj-2x2 


(P) 


-X1+X2  — ® 
•VX2  ^ -J 

xl*x2  - 0 


V 


-x^+x2  > 0 

xrx2  ^ 1 

x! *x2  > 0 


and  its  dual 

max  w 


<D) 


+2ul 


’“I’  U2 
1 1 
V u2 


< 0 
< -1 
< -2 


-u2  < 0 

•»S  i -1 
2 2., 
uj"u2  — ”Z 


u*  > 0,  i = 1,2;  k = 1,2 

The  primal  problem  (P)  has  an  optimal  solution  x = (0,2),  with  cx  = -4 
whereas  the  dual  problem'  (D)  is  infeasible.  This  is  due  to  the  fact  that 
q\q**  = [2}  and  X 2 = 0,  U2  = 0,  i.e.,  the  regularity  condition  is  violated. 


x e R, 


— xj+5C2  - 0 


xl"x2  - 1 


* U2 


u e R, 


2 2 . 
-V°2  < -1 

2 2 

V“2  * -2 


Here 


i 


5.  The  Convex  Hull  of  a Disjunctive  Set 
Having  described  the  family  of  all  valid  Inequalities,  one  Is  of 
course  Interested  In  identifying  the  strongest  ones  among  the  latter,  l.e., 
the  facets  of  the  convex  hull  of  feasible  points  of  a disjunctive  program. 
If  we  denote  the  feasible  set  of  a DP  by 


„n  V 


h h \ 

Ax2*o  U . 
*>o  h 


then  for  a given  scalar  the  family  of  inequalities  ax  > satisifed 
by  all  xeF,  l.e.,  the  family  of  valid  inequalities  for  the  DP,  is  obviously 
isomorphic  to  the  family  of  vectors  a € F^  y where 

F(a  ) “ (yeRn  I y*  > oQ>  VxcF}  , 

In  the  sense  that  ax  > a.  is  a valid  inequality  if  and  only  if  a c F , .. 

- 0 (aQ) 

In  view  of  its  relationship  with  ordinary  polar  sets,  we  call 
F^!  , the  reverse  polar  of  F.  Indeed,  the  ordinary  polar  set  of  F is 

(<v 

F°  = {y  e Rn  | yx  < 1,  Vx  e F}  , 

and  if  we  denote  by  F°^  ^ the  polar  of  F scaled  by  aQ,  (i.e.,  the  set 

0 ^ © 

obtained  by  replacing  1 with  a„  in  F),  then  F,  • -F,  ... 

0 (aQ)  (-aQ) 

The  size  (as  opposed  to  the  sign)  of  a^  is  of  no  interest  here.  There- 
fore we  will  distinguish  only  between  the  3 cases  a^  > 0 (or  Oq  = 1), 

= 0 and  Oq  < 0 (or  aQ  = -1).  (When  the  sign  of  ar^  makes  no  difference 
or  is  clear  from  the  context,  we  will  simply  write  F^ .)  For  a q < 0,  as 
mentioned  above,  F^  is  (the  negative  of)  an  ordinary  polar  set,  whose 
properties  are  described  in  the  literature.  The  most  interesting  case  for 
us,  however,  is  oq  > 1,  since  this  is  the  only  case  when  the  inequality 
ffx  > aQ  cuts  off  the  point  x = 0.  This  is  why  we  need  the  concept  of 


reverse  polars. 
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For  an  arbitrary  set  StRn  we  will  denote  by  cl  S,  conv  S,  cone  S, 
int  S and  dim  S,  the  closure,  the  convex  hull,  the  conical  hull,  the  interior 
and  the  dimension  of  S,  respectively.  For  a polyhedral  set  S£Rn  we  will 
denote  by  vert  S and  dir  S the  set  of  vertices  (extreme  points)  and  the  set 
of  extreme  direction  vectors  of  S,  respectively.  For  definitions  and  back- 
ground material  on  these  and  related  concepts  (including  ordinary  polar  sets), 
the  reader  is  referred  to  [ 37 ] or  [38]  (see  also  [28]). 

In  [10]  we  showed  that  while  some  of  the  basic  properties  of  polar  sets 
carry  over  to  reverse  polars,  others  can  only  be  recovered  in  a modified  form. 
In  the  first  category  we  mention  (a)  (XS)^  = ^ S^;  (b)  S CT  =»  s/; 

(c)  (S  U t/  = S^DT^,  properties  which  follow  from  the  definitions.  In  the 
second  one  we  state  a few  theorems,  which  are  from  [10]  (see  also  f 1 1 > . 

Theorem  5.1.  (i)  If  Oq  < 0,  then  j4  0 and 

0 e int  cl  conv  S <=**=>  is  bounded 


(ii)  If  a0  > 0,  then 


0 e cl  conv  S <=>  Sw  **  0 <=>  Sv  is  bounded. 


Proof,  (i)  follows  from  the  corresponding  property  of  the  ordinary 

nolar  S°  of  S and  the  fact  that  S^t  . = -S^  v. 

(a0)  («0> 

(ii)  For  oig  > 0,  if  S*  i4  0 there  exists  > e Rn  such  that  xy  > 


Vx  € cl  conv  S.  But  0 • y < o^,  hence  0 cl  conv  S.  Thus,  if  0 e cl  conv  S, 

then  S ^ = 0;  and  hence  is  bounded.  Conversely,  if  0 i cl  conv  S,  there 

exists  a hyperplane  ax  = a ^ separating  0 from  cl  conv  S,  i.e.,  such  that  (oq>0  and) 

ax  > Oq,  Vx  e cl  conv  S;  which  implies  aeS^,  i.e.,  / 0.  It  also  implies 

JL  JL 

XaeS  , VX>1,  i.e.,  S1f  is  unbounded.  |] 

From  this  point  on,  we  restrict  our  attention  to  sets  S whose  convex  hull  is 
polyhedral  and  pointed.  For  the  disjunctive  set  F,  this  condition  is  satisfied 
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if  Q is  finite.  Most  of  the  results  carry  over  to  the  general  case,  but 
proofs  arc  simpler  with  the  above  assumptions. 

Theorem  r>.2.  If  cl  conv  S is  polyhedral,  so  is  S^. 

Proof . Let  u^,...,u^  and  v^,...,v  be  the  vertices  and  extreme  direction 

vectors,  respectively,  of  cl  conv  S.  Then  for  every  yeS  there  exist  scalars 

P 

^ 2:  °>  i = P-  > 0,  j = l,...,q,  with  £ \i  = 1 , such  that 

i=l 

P q 

y = Eu.X.  + Ev.u., 

i-1  1 1 .1-1  J ’ 

and  it  can  easily  be  seen  that  for  arbitrary  x,  xy  > c*q,  V yeS,  if  and  only 
if  xu.  a(),  i - l,...,p,  and  xv^  >0,  j = l,...,q.  Thus 


x e R 


n|XUi  — G0  > 1 = 1 » • • • » P ] 


lxvj  > 0 > j = 1,. . .,qj 


i.e.,  S is  poTyhedral.ll 

y 

The  next  result  describes  the  crucial  involutory  property  of  polar  and 
reverse  polar  sets. 

Theorem  r>.').  Assume  / 0.  'Jhen 


cl  conv  S + cl  cone  S 


cl  cone  S 


cl  conv  (S  U [0]) 


if  c*q  > 0 

if  = ° 
if  a,.  < 0. 


Proof.  S#//  = {x  e Rn  | xy  > Vy  e S^J 


xy  > ctq  for  all  y s.t. 
x eRn  u y > <*0»  i * 1,...,P 


v±y  > °»  f * i>--..q 


1 


But  xy  > c<Q  is  a consequence  of  the  system  u^y  > Oq,  i ■=  and 

vAy  > 0,  i = (consistent,  since  S#  + (9),  if  and  only  if  there  exists 

a set  of  0^  > 0,  i = 1 , . . . ,p , > 0,  i = !,•••, q,  such  that 


(5.1)' 


x = E 0.  u,  + E a.v  , 
i=l  11  i=l 


.L*la0  ^ a0 

i=l 


Since  is  polyhedral,  so  is  S ^ . Thus  is  the  closed  set  of  points 

x e Rn  of  the  form  (5.1)  with  9^  > 0,  i = l,...,p  > 0,  i = l,...,q,  and 

' > 1 if  oQ  > 0 
P 

E 0.  = ( > 0 if  o:n  = 0 

i=l  1 U 

<1  if  a.  < 0. 

But  these  are  precisely  the  expressions  for  the  three  sets  claimed  in 
the  theorem  to  be  equal  to  in  the  respective  cases.  || 

Corollary  5.3.1.  cl  conv  S = Ds^^. 

Proof.  Follows  immediately  from  the  proof  of  Theorem  5.3,  where 

x e S,..  flS . ..  corresponds  to  £ 9 = l.|| 

\ * / \ *•)  ^ *■ 


Example  5.1.  Consider  the  disjunctive  set 


"X1  “x2  - “2 


,2  X1 


x2  - ® 


illustrated  in  Fig.  5.1(a) 


the  sets 


Xj>1  V x^>l 


Its  reverse  polars  for  q-q  = 1 and  «= 
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shown  in  Fig.  5.1(b)  and  (c). 

Finally,  the  sets  F##  corresponding  to  Qq  - 1 and  <*  = -1  (shown  in 

Fig.  5.2(a),  (b))  are 


X2  — 0 J 

and  their  intersection  (shown  in  Fig.  5.2(c))  is 


The  next  theorem  is  needed  to  prove  some  other  essential  properties  of 
reverse  polars. 

Theorem  5.4.  S = S . 

Proof.  If  cTq  < 0,  this  follows  from  the  corresponding  property  of 
ordinary  polars.  If  c*q  > 0 and  0 € cl  conv  S,  then  = 0,  = Rn,  and 

Sm  = 0 = S#.  Finally,  if  > 0 and  0 l cl  conv  S,  then 

Sm  = cl  (conv  S + cone  S)^  (from  Theorem  5.3) 

= {y  6 Rn  | xy  > Oq»  Vx  e cl  (conv  S + cone  S ) } 

= {y  e Rn  | xy  > ar0,  Vx  e S}  = S#.  || 
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The  above  results  can  be  used  to  characterize  the  facets  of  cl  conv  S. 

To  simplify  the  exposition,  we  will  assume  that  S is  a full-dimensional 
pointed  polyhedron,  which  implies  that  is  also  full-dimensional.  For  the 

general  case  the  reader  is  referred  to  [10]. 

We  recall  that  an  inequality  rrx  > tt^  defines  a facet  of  a convex  set  C 
if  ttx  > t^,  VxeC,  and  nx  = tt^  for  exactly  d affinely  independent  points  x 
of  C,  where  d = dim  C.  The  facet  of  C defined  by  nx  > nQ  is  then 
{x e C | nx  = fig];  but  as  customary  in  the  literature,  for  the  sake  of  brevity 
we  call  nx  > n^  itself  a facet. 

We  proceed  in  two  steps,  the  first  of  which  concerns  the  facets  of  S^. 


Theorem  5.5.  Let  dim  S = n,  and  t4  0.  Then  ax  > Bq  is  a facet  of 
S##  if  and  only  if  a ^ 0 is  a vertex  of  . 

Proof.  From  Theorem  5.4,  a e Rn  is  a vertex  of  if  and  only  if 


f 

uy  > a„>  Vu  evert  j 

= < 

. yt»" 

vy  > 0 , Vv  e dir  j 

and  a satisfies  with  equality  a subset  of  rank  n of  the  system  defining  S^. 

Further,  a ^ 0 if  and  only  if  this  subset  of  inequalities  is  not  homogeneous 

(i.e.,  at  least  one  right  hand  side  coefficient  is  aQ  t 0). 

On  the  other  hand,  ax  > aQ  is  a facet  of  if  and  only  if  (i)  Qx  > a0> 

Vx  e , i.e.,  aeS^;  and  (ii)  ax  = aQ  for  exactly  n affinely  independent 

$4 

points  of  S'  . But  (ii)  holds  if  and  only  if  a u = Q'q  for  r vertices  u of 
S*#,  and  av  = 0 for  s extreme  direction  vectors  v of  S^,  with  r > 1 
(since  a / 0)  and  r + s > n,  such  that  the  system  of  these  equations  is  of 


J 


rank  n. 


Thus  the  two  sets  of  conditions  (for  ax  > Og  to  be  a facet  of  and 
for  a to  be  a vertex  of  S^)  are  identical.  || 

By  arguments  similar  to  the  above  proof,  one  shows  that  ox  > 0 is  a 

JLJL 

facet  of  S™  if  and  only  if  a is  an  extreme  direction  vector  of  S . Unlike 
for  Oq  Vs  0,  the  homogeneous  inequality  ax  > 0 is  a facet  of  if  and  only 

if  it  is  also  a facet  of  [of  S^],  due  to  the  fact  that  every  extreme 

direction  vector  of  is  also  an  extreme  direction  vector  of 

[of  S^q) ] , and  vice-versa. 

Theorem  5.6.  Let  dim  S = n,  and  a^  / 0.  Then  ax  > a^  is  a facet  of 


cl  conv  S if  and  only  if  it  is  a facet  of  S 


## 


Proof,  (i)  If  aQ  > 0,  the  half space  ax  > aQ  contains  cl  conv  S if  and 

only  if  it  contains  cl(conv  S + cone  S) . If  ag  < 0,  the  halfspace  ax  > ap 

contains  cl  conv  S if  and  only  if  it  contains  cl  conv  (SU{0)>.  From  Theorem  4.3, 

in  both  cases  ax  > Oq  is  a supporting  halfspace  for  cl  conv  S if  and  only 

if  it  is  a supporting  halfspace  for  . 

(a0) 

(ii)  Next  we  show  that 


{x  e cl  conv  S | a x = ag}  = {x  e ^ |ax  = aQ}  . 

The  relation  £ follows  from  cl  conv  S £ . (Theorem  5.3).  Tc  show  the 

(q0) 
u ## 

converse,  assume  it  to  be  false,  and  let  x e S.  Acl  conv  S satisfy  ax  = a«. 

(aQ)  0 

From  Theorem  5.3,  x = \u  for  some  u e cl  conv  S,  and  1 > 1 if  Oq  > 0, 

0 < X < 1 if  Oq  < 0.  In  each  case,  a x *=  ag  implies  au  * (1/X)  ax  < ag  for 

some  u e cl  conv  SCS^  . , contrary  to  the  assumption  that  a*  > a.^  x e v 

(a0)  0 
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By  an  argument  similar  to  the  above  proof  one  can  show  that  if  a x > 0 


for  both  ar 


1 and 


1 , 

* [ 

| x e R 

heQ  ^ 

is  a facet  of  cl  conv  S,  then  it  is  a facet  of  S ^ ^ Uq 

erg  = -1.  The  converse,  however,  is  not  true,  i.e.,  a x > 0 can  be  a facet 
of  both  and  S ^_iy  without  being  a facet  of  OS ^iy 

We  are  now  ready  to  characterize  the  facets  of  the  convex  hull  of  the 
disjunctive  set 

/ 

i . h h 

[A  x >aQ 

x > 0 

where  Q is  assumed  to  be  finite,  and  F to  be  full-dimensional  (for  the 
general  case  see  [10]). 

Theorem  5.7.  ax  > oQ,  with  M,  is  a facet  of  cl  conv  F if  and 
only  if  a i 0 is  a vertex  of  the  polyhedron 


. h li 

where  Q is  the  set  of  those  h e Q for  which  the  system  A x > aQ,  x > 0,  is 

consistent. 

Proof.  From  Theorem  3.1,  the  set  F^  ^ = {y  e Rn  | xy  > Oq,  Vx  sF}  is  of 
the  form  claimed  above.  The  rest  is  a direct  application  of  Theorems  5.5 
and  5.6.  II 


\ 

h h * 

y > 6 A , hcQ 

yeRn 

for  some  0 > 0,  h eQ  ) 

satisfying  0hag  > aQ 

/ 

As  for  the  case  qiq  = 0,  from  the  above  conxnents  it  follows  that  if 
a x > 0 is  a facet  of  cl  conv  F,  then  a t 0 is  an  extreme  direction  vector 
of  F^  \ for  all  Oq»  The  converse  is  not  true,  but  if  a ^ 0 is  an  extreme 


(cro) 
direction 


vector  of  Y^.  ^ for  some  qq  (hence  for  all  Og),  then  or  x > 0 is 

0 11 
either  a facet  of  cl  conv  F,  or  the  intersection  of  two  facets,  or  x > ag  and 

a2x  > Og,  with  o1  + a2  = a and  Og  = -a2  t 0 (see  [10]  for  the  details). 
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Since  the  facets  of  cl  conv  F are  vertices  (in  the  nonhomogeneous  case) 
or  extreme  direction  vectors  (in  the  homogeneous  case)  of  the  convex  poly- 
hedron  F , they  can  be  found  by  maximizing  or  minimizing  some  properly  chosen 
linear  function  on  , i.e.,  by  solving  a linear  program  of  the  form 

min  gy 


P1<8’V 

or  its  dual 


y - 0hAh  > 0 

fih  h ^ 

6 a0  ^ ao 

0h  > 0 


h e 


* 

Q 


max 


* 


?2 ^ »®Q ) 


E ^ § 8 

heQ 


From  Theorem  5.1, 
feasible;  whereas  if 


h_h  .h_h  - _ , * 

a0§Q  “ A S < 0 heQ 

> °»  > o 

if  CfQ  < 0 then  F^  ) J4  0,  i.e.,  P*(g,aQ)  is 
* 0 

> 0,  then  P^(g,o,Q)  is  feasible  if  and  only 


always 

if 


0 i cl  conv  F.  This  latter  condition  expresses  the  obvious  fact  that  an 

inequality  which  cuts  off  the  origin  can  only  be  derived  from  a disjunction 

which  itself  cuts  off  the  origin. 

Two  problems  arise  in  connection  with  the  use  of  the  above  linear 

programs  to  generate  facets  of  cl  conv  F.  The  first  one  is  that  sometimes 

only  Q is  known,  but  Q is  not.  This  can  be  taken  care  of  by  working  with 

Q rather  than  Q . Let  P.  (g,o  ) denote  the  problem  obtained  by  replacing  Q 

K 0 

ic 

with  Q in  P^gjOg),  k = 1,2.  It  was  shown  in  [10],  that  if  ?2^>aQ^  has  an 
optimal  solution  | such  that 


(1q  - °,  th  «*  0)  =»  h c Q* 


j\J (Xv*\t  j 


F 


jt  r\ 


/VCX/v  \Jc 


Qo^ 


J' 


3. GIVEN  THE  OPTION  OF  PURCHASING  ANY  ONE  OF  TWO  STOCKS  A AND  B 
WITH  THE  FOLLOWING  RISK-RETURN  CHARACTERISTICS.  JOHN  SAYS 
HE  WILL  BUY  STOCK  A,  FROM  THIS  DATA  I CAN  CONCLUDE  THAT: 


EXPECTED  RETURN 

RISK 

STOCK 

A 

57. 

87. 

STOCK 

B 

67. 

9 7. 

A. 

JOHN 

IS 

A RISK  LOVER 

B. 

JOHN 

IS 

RISK  ADVERSE. 

C. 

JOHN 

IS 

NEUTRAL  TOWARDS  RISK 

D. 

JOHN 

IS 

IRRATIONAL 

E. 

NONE 

OF 

THE  ABOVE. 

1 


4. NOW  SUPPOSE  THAT  JOHN  HAS  THE  CHOICE  OF  ANY  ONE  OF  THREE 
STOCKS  A » B AND  C.  A AND  B HAVE  THE  SAME  RISK-RETURN 
CHARACTERISTICS  AS  IN  QUESTION  3.  STOCK  C HAS 

EXPECTED  RETURN  OF  5%  AND  RISK  OF  9%.  WHICH  STOCK  WILL  JOHN  CHOOSE 


A. 

B. 

C. 

D. 

E. 


STOCK  A 


STOCK  B 
STOCK  C 

UNDETERMINED  WITH  GIVEN  INFORMATION 

HE  WOULD  BE  INDIFFERENT  BETWEEN  STOCKS  B AND  C 


then  every  optimal  solution  of  P^Cg,^)  is  an  optimal  solution  of  P^g,^). 
Thus,  one  can  start  by  solving  P2(g,ag).  If  the  above  condition  is  violated 

lAr 

for  some  h «Q\Q  , then  h can  be  removed  from  Q and  P2(g »cq)  solved  for  Q 
redefined  in  this  way.  When  necessary,  this  procedure  can  be  repeated. 

The  second  problem  is  that,  since  the  facets  of  cl  conv  F of  primary 
interest  are  the  nonhomogeneous  ones  (in  particular  those  with  > 0,  since 
they  cut  off  the  origin),  one  would  like  to  identify  the  class  of  vectors  g 
for  which  P^(g,ag)  has  a finite  minimum.  It  was  shown  in  [10],  that  P^(g,<jQ) 
has  a finite  minimum  if  and  only  if  \g  eel  conv  F for  some  \ > 0;  and  that, 
should  g satisfy  this  condition,  ax  > Oq  is  a facet  of  cl  conv  F (where  F 
is  again  assumed  full-dimensional)  if  and  only  if  a = y for  every  optimal 
solution  (y,§)  to  P1(g,a0). 

As  a result  of  these  characterizations,  facets  of  the  convex  hull  of 
the  feasible  set  F can  be  computed  by  solving  the  linear  program  P^(g,ag)  or 
its  dual.  If  the  disjunction  defining  F has  many  terms,  like  in  the  case 
where  F comes  from  the  disjunctive  programming  formulation  of  a 0-1  program 
with  a sizeable  number  of  0-1  conditions,  P^(g,aQ)  is  too  large  to  be  worth 
solving.  If,  however,  F is  made  to  correspond  to  a relaxation  of  the  original 
zero-one  program,  involving  zero-one  conditions  for  only  a few  well  chosen 
variables,  then  Pj(g»org)  or  its  dual  is  practically  solvable  and  provides  the 
strongest  possible  cuts  obtainable  from  those  particular  zero-one  conditions. 

On  the  other  hand,  since  the  constraint  set  of  P2(g,Og)  consists  of  |q  | 
more  or  less  loosely  connected  subsystems,  one  is  tempted  to  try  to  approximate 
an  optimal  solution  to  P^g.Og)  — and  thereby  to  P1(g,OQ)  — by  solving  the 
subsystems  independently.  Early  computational  experience  indicates  that 


these  approximations  are  quite  good. 

We  now  give  a numerical  example  for  a facet  calculation 
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Example  5.2.  Find  all  those  facets  of  cl  conv  F which  cut  off  the 

2 

origin  (i.e.,  all  facets  of  the  form  or  x > 1),  where  FCR  is  the  disjunctive 
set 

F - Tl  V F2  V F3  V F4  , 

with 

Fx  = {x  | -Xj^  + 2x2  > 6,  0 < Xj  < 1,  x2  > 0} 

F2  = [x  | 4xx  + 2x2  >11,  1 < xx  <2.5,  x2  > 0) 

F3  - (x  | -xx  + x2  > -2,  2.5  < x1  < 4,  x2  > 0] 

F4  “ (x  I x!  + x2  ^ 6»  4 5 XL  < 6»  x2  ^ °} 

(see  Fig.  5.3). 

After  removing  some  redundant  constraints,  F can  be  restated  as  the 
2 

set  of  those  xtR+  satisfying 

| 4xj  + 2x?  > 11  | 

{-x.  + 2x2  > 6}  V / ) V {x  + x > 6} 

rxi+  x2  — ~2  j 

and  the  corresponding  problem  P^(g,l)  is 


min  g]yL  + g2y2 


+ 0 


y2  - 2e; 

y2 

y2 

6el 


- 402  + 02 

2 2 
-261-  02 


1102  - 202 


0 

0 

0 

0 

0 

0 

1 

1 

1 


12  2 3 

®1»  ®1»  02*  01  ^ ° 


32 


Solving  this  linear  program  for  g - (1,1),  yields  the  optimal  points 

(y I®)  * > 6**^9* 9*^)  1 and  (y»0)  **  ^3* *3*  * 6^ **9 **9 * which  have  the  same 

y-component:  (^,^).  These  points  are  optimal  (and  the  associated  y is 

unique)  for  all  g > 0 such  that  gx  < 5g2.  For  gL  *>  5g2>  in  addition  to  the 

above  points,  which  are  still  optimal,  the  points  (y;9)  -(t,t 

\o  6 o 9 18  6/ 

and  (y;0)  = which  again  have  the  same  y-component 

y = also  become  optimal;  and  they  are  the  only  optimal  solutions  for 

all  g > 0 such  that  g^  > 5g2« 

We  have  thus  found  that  the  convex  hull  of  F has  two  facets  which  cut 
off  the  origin,  corresponding  to  the  two  vertices  y*  = 0,^)  and  y^  = (^, ~J 
of  F#(l): 

1 ^ 1 

3X1+3X2  ^ 1 
6 X1  + 6 x2  ^ 1 


is  + 


6 Xi  f ( X.Z  - 1 


(0,^  If  ,0)  (i^O)  (*,ol 


Fig.  5.3 
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6.  Facial  Disjunctive  Programs 

In  this  section  we  discuss  the  following  problem  [10].  Given  a disjunctive 
program  in  the  conjunctive  normal  form  (2.2),  is  it  possible  to  generate  the 
convex  hull  of  feasible  points  by  imposing  the  disjunctions  j cS  one  by  one, 
at  each  step  calculating  a "partial"  convex  hull,  i.e.,  the  convex  hull  of  the 
set  defined  by  the  inequalities  generated  earlier,  plus  one  of  the  disjunctions? 

For  instance,  in  the  case  of  an  integer  program,  is  it  possible  to 
generate  the  convex  hull  of  feasible  points  by  first  producing  all  the  facets 
of  the  convex  hull  of  points  satisfying  the  linear  inequalities,  plus  the 
integrality  condition  on,  say,  x^;  then  adding  all  these  facet-inequalities 
to  the  constraint  set  and  generating  the  facets  of  the  convex  hull  of  points 
satisfying  this  amended  set  of  inequalities,  plus  the  integrality  condition  on 
X2;  etc.  Hie  question  has  obvious  practical  importance,  since  calculating 
facets  of  the  convex  hull  of  points  satisfying  one  disjunction  is  a consider- 
ably easier  task,  as  shown  in  the  previous  section,  than  calculating  facets 
of  the  convex  hull  of  the  full  disjunctive  set. 

The  answer  to  the  above  question  is  negative  in  general,  but  positive 
for  a very  important  class  of  disjunctive  programs,  which  we  term  facial. 

The  class  includes  (pure  or  mixed)  0-1  programs. 

The  fact  that  in  the  general  case  the  above  procedure  does  not  produce 
the  convex  hull  of  the  feasible  points  can  be  illustrated  on  the  following 
2-variable  problem. 

Example  6.1.  Given  the  set 
2 

F0  = (x  e R | -2x^  + 2-<2  < 1,  2x1  - 2x2  <1,  0 < x1  < 2,  °<*2—  2J 


find  F ■ cl  conv  (Fq  0 {x  | x^,X2  integer]) 


Denoting 


F^  = cl  conv  (Fq  0 {x  | x^  integer}),  F 2 ° cl  conv  (i?in{x|x2  integer}), 

the  question  is  whether  F2  ' F.  As  shown  in  Fig.  6.1,  the  answer  is  no, 
since 


/ 


2x 

2x. 

■2x 

2x 


+ 3x2  > 0 
+ x2  > -2 
- 3x2  > ~2j 


while 


F 


/ 

\ 

-x,  + x„  = 0 

1 2 

< X 

° < X!  < 2 ) 

0 < x„  <2 

k 

— 2 — 

If  the  order  in  which  the  integrality  constraints  are  imposed  is 


reversed,  the  outcome  remains  the  same. 


t - JlcxK*  (tni* I*Z iJp+l)  ^ ^ V 
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Cons  ider  a disjunctive  program  stated  in  the  conjunctive  normal  form 
(2.2),  and  denote 

Fq  - {x  cRn  | Ax  > aQ,  x > 0}  . 

The  disjunctive  program  (and  its  constraint  set)  is  called  facial  if 
every  inequality  d^x  > d^  that  appears  in  a disjunction  of  (2.2),  defines 
a face  of  FQ;  i.e.,  if  for  all  i eQj,  j e S,  the  set 

FQn{x  | d*x  > d?Q} 

is  a face  of  Fq.  (A  face  of  a polyhedron  P is  the  intersection  of  P with 
some  of  its  boundary  planes.) 

The  class  of  disjunctive  programs  that  have  the  facial  property  includes 
the  most  important  cases  of  disjunctive  programming,  like  the  0-1  programming 
(pure  or  mixed),  nonconvex  quadratic  programming,  separable  programming,  the 
linear  complementarity  problem,  etc.;  but  not  the  general  integer  programming 
problem,  as  illustrated  above.  In  all  the  above  mentioned  cases  the 
inequalities  d*x  > d ^ of  each  disjunction  actually  define  facets,  i.e., 

(d  - 1) -dimensional  faces  of  Fq,  where  d is  the  dimension  of  Fq. 

Another  property  that  we  need  is  the  boundedness  of  Fq.  Since  this  can 
always  be  achieved,  if  necessary,  by  regularizing  Fq,  its  assumption  does 
not  represent  a practical  limitation. 

Theorem  6.1.  Let  the  constraint  set  of  a DP  be  represented  as 

F - {*«F0|  V w‘x>d  ).  jes} 

UQj 

where 

Fq  = {x  e Rn  | Ax  > aQ , x > 0}  . 

For  an  arbitrary  ordering  of  S,  define  recursively 
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Fj  » conv 


[ v (p  n(*  I d1^  > ’i0')J  . i - i..  . ;* I . 

4 J 


i*Q 


If  F is  facial  and  Fq  is  bounded,  cf e \ 


F Is  I " conv  r * 


The  proof  of  this  theorem  (10’  uses  the  following  a. »i'  y re  >1 


Lemma  6.1.1.  Let  P^,...,Pr  bi  a finite  set  ^olytopts  ^bounned 
r 

polyhedra),  and  P = UP.. 

h=l 

Let  H+  = {x  c Rn  | dSc  < d^g}  be  an  *rh  naltspa  - 

H ■ {xcRn  | dix  * diJ}  its  defining  hyperpltne. 

If  P c H+,  then 

HOconv  P » -onv  ('d.i0)  . 


Proof.  Let  HO  conv  P j t (otherwise  the  Lenina  holds  trivially). 
Clearly,  (HOP)  £ (HD conv  P),  and  therefore 

conv  (HOP)  c conv  (Hflconv  P)  ■ Hfl  conv  P . 

To  prove  2,  let  u^,...,Up  be  the  vertices  of  all  the  poly  topes 

p , h = l,...,r.  Obviously,  p is  finite,  conv  P is  closed,  and 
n 

r 

vert  conv  P £ ( U vert  P,  ).  Then 
h°l 


xiHfl conv  P =» 


dlx 


iO 


and  x 


k=l 


with 

P 

E = 1,  Xfc  > °»  k = 1.....P  • 

k=l 

Further,  PcH+  implies  d^  < diQ,  k «*  1 We  claim  that  in  the  above 

expression  for  x,  if  Xk  > 0 then  diuk  = diQ.  To  show  this,  we  assume  there 
exist  Xk  > 0 ouch  that  diufe  < diQ. 


Then 
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i i p P 

d1*  - dx(  Z X.iO  < d Q(  Z X.) 

k-1  * K 1U  k-1  K 

" diO  * 

u 

a contradiction.  Hence  is  the  convex  combination  of  points  u cHHP,  or 
x « conv  (HOP).  || 


Another  relation  to  be  used  in  the  proof  of  the  theorem,  is  the  fact 
that  for  ai c-i  trary  S^,  CRn, 

- 1 ' -onv  (conv  S^Uconv  S2)  “ conv  (S^US2)  . 

Proof  01  Theorem  6.1.  For  j ■*  1,  the  statement  is  true  from  the 
definitions  and  the  obvious  relation 

V (F0fi[x  | d*x  > d 0))  - {xcF0|  V (dlx>d  )}. 

ieQl  ieQ! 

To  prove  the  theorem  by  induction  on  j,  suppose  the  statement  is  true 
for  j ■ l,...,k.  Then 

* conv  [ V (Fkn{x  | d*x  > d^})]  (by  definition) 

ieQkfl 

= conv  [ V ({x  | d£x  > d 0}  n conv  [x  cFQ  | V (d*x  > d,_),  J « l,...,k)) 
i€Qkfl  ieQj 

(from  the  assumption) 


= conv  [ V conv  ({x  | d*x  > d } 0 {x  eFQ  | V (d*x  > d ) , j - 1 k} ) ] 

ieQkfl  icQj 

(from  Lemma  6.1.1) 


" conv  [(  V (x  | d£x  > diQ})  0 {x  *F0  | V (d*x  > d Q)f  J - 1 k) ] 

ieQlcH  ieQj 

(from  (6.1)) 


= conv  Cx  c F Q | V (dix  > d Q),  j - 1 V*l}  , 

i£QJ 
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i.e.,  the  statement  is  also  true  for  j - k + 1.  || 

Theorem  6.1  implies  that  for  a bounded  facial  disjunctive  program  with 
feasible  set  F,  the  convex  hull  of  F can  be  generated  in  |s|  stages,  (where 
S is  as  in  (2.2)),  by  generating  at  each  stage  a "partial"  convex  hull, 
namely  the  convex  hull  of  a disjunctive  program  with  only  one  disjunction. 

In  terms  of  a 0-1  program,  for  instance,  the  above  result  means  that 
the  problem 

min  [ex  | Ax  > b,  0 < x < e,  x^  = Oor  1,  j ” 1,...  ,n}  , 

where  e = (1,...,1),  is  equivalent  to  (has  the  same  convex  hull  of  its 
feasible  points,  as) 

(6.2)  min  [ex  | Ax  > b,  0 < x < e,  a1*  > c*10,  UHj,  Xj  = 0 or  1,  j = 2,...,n} 

where  a^x  > i are  the  facets  of 

F^  ■ conv  {x  | Ax  > b,  0 < x < e,  - 0 or  1}  . 

In  other  words,  x^  is  guaranteed  to  be  integer-valued  in  a solution  of  (6.2) 
although  the  condition  ■ 0 or  1 is  not  present  among  the  constraints  of 

(6.2) .  A 0-1  program  in  n variables  can  thus  be  replaced  by  one  in  n-1 
variables  at  the  cost  of  introducing  new  linear  inequalities.  The  inequalities 
themselves  are  not  expensive  to  generate,  since  the  disjunction  that  gives 
rise  to  them  (x^  = 0 V x^  = 1)  has  only  two  terms.  The  difficulty  lies  rather 
in  the  number  of  facets  that  one  would  have  to  generate,  were  one  to  use  this 
approach  for  solving  0-1  programs.  However,  by  using  some  information  as  to 
which  inequalities  (facets  of  a "partial"  convex  hull)  are  likely  to  be 
binding  at  the  optimum,  one  might  be  able  to  make  the  above  approach  efficient 
by  generating  only  a few  facets  of  the  "partial"  convex  hull  at  each  iteration. 
This  question  requires  further  investigation.  For  additional  results  on 
facial  disjunctive  programs  see  [33],  [34], 
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7.  Dislunctive  Programs  with 


Explicit  Integrality  Constraints 


Hie  theory  reviewed  in  the  previous  sections  derives  cutting  planes 
from  disjunctions.  In  this  context,  0-1  or  integrality  conditions  are  viewed 
as  disjunctions,  and  the  disjunction  to  be  used  for  deriving  a cut  usually 
applies  to  the  basic  variables. 

In  this  section  we  discuss  a principle  for  strengthening  cutting  planes 
derived  from  disjunctions  in  the  case  when,  besides  the  disjunction  which 
applies  to  the  basic  variables,  there  are  also  integrality  constraints  on 
some  of  the  nonbasic  variables.  In  [14]  we  first  proved  this  principle  for 
arbitrary  cuts,  by  using  subadditive  functions,  then  applied  it  to  cuts 
from  disjunctions.  Here  we  prove  the  principle  directly  for  the  latter  case, 
without  recourse  to  concepts  outside  the  framework  of  disjunctive  programning. 

Let  a DP  be  stated  in  the  disjunctive  normal  form  (2.1),  and  assume 
in  addition  that  some  components  of  x are  integer-constrained.  In  order  for 
the  principle  that  we  are  going  to  discuss  to  be  applicable,  it  is  necessary 
that  each  A^x,  heQ,  has  a lower  bound,  say  bjj.  With  these  additional 
features,  and  denoting  by  J the  index  set  for  the  components  of  x (|j|  = n), 
the  constraint  set  of  the  DP  can  be  stated  as 


(7.1) 


A*x  > bj , i e Q 
x > 0 , 


(7.2) 


V (A*x  > aj) 

ieQ 


(7.3) 


where 


(7.4) 


Xj  integer,  jc^cj  , 


i.Q 
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Let 


Q = and  let  a*  stand  for  the  j-th  column  of  A1,  j « J, 


1 «Q. 


Theorem  7.1.  [14]  Define 

(7.5)  M - [m  eRq  | E m1  > 0,  mt  integer,  i e q]  . 

i«Q 

Then  every  x c Rn  that  satisfies  (7.1),  (7.2),  (7.3),  also  satisfies  the 
inequality 


(7.6) 

where 


(7.7)  aj 


E £ “0  » 
JeJ  J J 


' inf  max  0*[a*  + m.(aQ-b^)]  , J € 
meM  ieQ  L 

J max  0*a*  Jej\j^  * J2 

V i*Q  J 


and 

(7.8)  an  ■ min  0ia^  . 

0 i«Q 


To  prove  this  theorem  we  will  use  the  following  auxiliary  result. 

Lemma  7.1.  Let  m^  c M,  m^  *=  (m^),  j e Jj.  Then  for  every  x e Rn 
satisfying  (7.3),  either 

(7.9)  E miJXj  “ ° * V1  eQ 

J*Jj 

or 

(7.10)  V ( E m x > 1) 

ieQ  JeJ1  11  3 i, 


holds . 
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Proof.  If  the  statement  is  false,  there  exists  x satisfying  (7.3)  and 
such  that 

E E < 0 * 

ieQ  J€JX  J 

On  the  other  hand,  from  x > 0 and  the  definition  of  M, 


E 2 m x > 0 , 
ieQ  JeJ1  J J 

a contradiction.  || 

Proof  of  Theorem  7.1.  We  first  show  that  every  x which  satisfies  (7.1), 
(7.2)  and  (7.3),.  also  satisfies 

<7-2,) 


for  any  set  of  cM,  j e J^.  To  see  this,  write  (6.2  ')  as 

(7.2")  V[E  a*x  + (aj.-bb  E m x > aj]  . 

ieQ  jeJ  J J u u jeJj  lJ  J 

From  Lemma  7.1,  either  (7.9)  or  (7.10)  holds  for  every  x satisfying  (7.3). 

If  (7.9)  holds,  then  (7’.2")  is  the  same  as  (7.2)  which  holds  by  assumption. 

If  (7.10)  holds,  there  exists  keQ  such  that  E 1\jxj  = 1 + X for  some  X > 0. 

jeJj 

But  then  the  k-th  term  of  (7.2")  becomes 


_ k ^ . k . k , k s 

E ajxj  - bo  ' x(ao-V 

JCJ 

k k 

which  is  satisfied  since  X(Sq  - b^)  > 0 and  x satisfies  (7.1).  This  proves 
that  every  feasible  x satisfies  (7.2/). 

Applying  to  ( 7 . 2 7 ) Theorem  3.1  then  produces  the  cut  (7.6)  with 


coefficients  defined  by  (7.7),  (7.8).  Taking  the  infinum  over  M is  justified 


i 


by  the  fact  that  (7.6)  is  valid  with  an  as  in  (7.8),  as  in  (7.7)  for 


j c J, , and 


a.  - max  61 [a*  + m (a*  - b*)] 
3 ieQ  J J 


for  j e Jp  for  arbitrary  e M.  || 

Corollary  7.1.1.  [14],  Let  the  vectors  a*,  ieQ,  satisfy 


(7.11) 


i,  i , i,  , ii\  . 
o (a0-bQ)  = 1,  aaQ  >0  . 


Then  every  x e Rn  that  satisfies  (7.1),  (7.2)  and  (7.3),  also  satisfies 


(7.6') 


where 


E > 1 » 


(7.7 ')  Pj  - < 


i i ^ 
o a + m 

min  max  

i i 

meM  ieQ  a Sq 


Oia^ 


max 


ieQ  a a 


i i * 


j e J. 


j * J, 


\ 

Proof.  Given  any  a1,  ieQ,  satisfying  (7.11),  if  we  apply  Theorem  7.1 
by  setting  01  = (oVc/a^) , ieQ,  in  (7.7)  and  (7.8),  we  obtain  the  cut  (7.6'), 
with  Pj  defined  by  (7.7'),  j e J.  || 

Note  that  the  cut-strengthening  procedure  of  Theorem  7.1  requires,  in 
order  to  be  applicable,  t.he  existence  of  lower  bounds  on  each  component  of 
A*x,  Vi  eQ.  This  is  a genuine  restriction,  but  one  that  is  satisfied  in  many 


practical  instances.  Thus,  if  x is  the  vector  of  nonbasic  variables  associated 
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r> 

with  a given  basis,  assuming  that  A^x  is  bounded  below  for  each  icQ 
amounts  to  assuming  that  the  basic  variables  are  bounded  below  and/or 
above.  In  the  case  of  a 0-1  program,  for  instance,  such  bounds  not  only 
exist  but  are  quite  tight. 

Example  7.1.  Consider  again  the  mixed-integer  program  of  example  3.1 
(taken  from  [35]),  and  assume  this  time  that  and  x^  are  0-1  variables 

rather  than  just  integer  constrained,  i.e.,  let  the  constraint  set  of  the 
problem  be  given  by 


X1 

* .2  + ,4(-x3)  + 

1*3 (-x^) 

- .01(-x5) 

+ .07(-x6) 

X2 

= .9  - .3(-x3)  + 

•4(-x4) 

- .04(-x5) 

+ .1  (-x6) 

> 0, 

j = 1,. . . ,6;  Xj  = 

: 0 or  1, 

j = 1,2; 

Xj  integer 

Biis  change  does  not  affect  the  Gomory  cuts  or  the  cuts  obtainable 
from  extreme  valid  inequalities  for  the  group  problem,  which  remain  the 
same  as  listed  in  example  3.1. 

Now  let  us  derive  a cut,  strengthened  by  the  above  procedure,  from 
the  disjunction 

x.  > (0 

> v {*2  > 1)  • 
x2  > °J 

Since  x^,  -X2  and  X2  are  bounded  below  by  0,  -1  and  0 respectively,  we 

have 


1 2 

Applying  Corollary  7.1.1,  we  choose  a = (4,1),  o = 1,  which  is  easily- 


seen  to  satisfy  (7.11).  Since  Q has  only  2 elements,  the  set  M of  (7,5), 


becomes 


M - {m  - (n^.n^)  | m1  + m2  > 0;  n^,  m2  integer} 
and,  since  at  the  optimum  we  may  assume  equality, 

M = [m  = (m^-n^)  | integer}  . 

The  coefficients  defined  by  (7 . 7 7 ) then  become 


min 
m^  integer 


max 


{ 


4 X (-.4)+ 1 x (-•3)  + m1  lx.3-mi 


4 x (-.2)+  1 X .9 


lx  - 1 / 


-7  (with  « 


min 
m^  integer 


max 


r 4 X (-1.3)+  1 x .4  + m^  1 X (-.4)  -m.^  N 
1 4 x (-.2)+  1 x .9  ’ 1 X .1  J 


-24  (with  m 


P5 

^6 


max 


max 


f 4 x .01+  1 X (-.04)  1 X .04  1 

l 4 x (-.2)  + 1 X .9  ’ 1 X .1  J 

v 


r 4 x (-.07)  + 1 x .1  1 X (-.1)  ) ^ , 

l 4 X (-.2)+ 1 X .9  * lx.l  J 


and  the  cut  is 


— 7x^  — ■ 24x^  + .4x,-  — x6  1 > 

which  has  a smaller  coefficient  for  x^  (and  hence  is  stronger)  than  the 
cut  derived  in  example  3,1. 

In  the  above  example,  the  integers  nu  were  chosen  by  inspection. 
Derivation  of  an  optimal  set  of  na  requires  the  solution  of  a special  type 
of  optimization  problem.  Two  efficient  algorithms  are  available  [14]  for 
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doing  this  when  the  multipliers  o are  fixed.  Overall  optimization  would 
of  course  require  the  simultaneous  choice  of  the  o‘  and  the  m^ , but  a good 
method  for  doing  that  is  not  yet  available. 

The  following  algorithm,  which  is  one  of  the  two  procedures  given  in 
[14],  works  with  fixed  o*,  i eQ.  It  first  finds  optimal  noninteger  values 
for  the  m^,  i cQ,  and  rounds  them  down  to  produce  an  initial  set  of  integer 
values.  The  optimal  integer  values,  and  the  corresponding  value  of  p , 
are  then  found  by  applying  an  iterative  step  k times,  where  k < |q | - 1, 

|Q | being  the  number  of  terms  in  the  disjunction  from  which  the  cut  is 
derived. 


Algorithm  for  calculating 


e J^,  of 


Denote 


(7.12) 


ori  = a , 


. , i i. -1 

Xi  = (o  ao> 


(7.13) 


E a.  / E r- 
icQ  ieQ  Ai 


Calculate 


(7.14) 


* X 

mi  = x ' «i  ' 


i eQ  • , 


“fc 

set  m.  = [m, ] , i eQ,  define  k = - £ [m.],  and  apply  k times  the  following 

1 icQ  1 


I terative  Step . Find 


X8(c*8  + ros  + !)  = ">in  X^o^  + 1^+1) 

icQ 


and  set 


m •-  m + 1 , 
8 s 


mi  *”  mi  * 


i eQ\{s] 
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This  algorithm  was  shown  in  [14]  to  find  an  optimal  set  of  m^  (and 

the  associated  value  of  g.)  in  k steps,  where  k = - E [m  ] < |Q  | - 1 . 

J ieQ 

Example  7.2.  Consider  the  integer  program  with  the  constraint  set 

X1  c 6 + F (“x5 ) “ (_x5 ) 'g  (_X7 ) 

x2  " | + £ + I (-x6)  " £ <"x7> 

3 2 4 1 

x3  - f " £ (-x5)  + | (-x6>  - £ <‘x7> 

x4  = £ + I ("x5)  + I (-x6)  " £ (_x7) 

X,  + X_  + X.  + X,  > 1 

1 2 3 4 - 

Xj  = 0 or  1,  j - 1,...,4;  x^  > 0 integer,  j - 5,6,7. 

We  wish  to  generate  a strengthened  cut  from  the  disjunction 
x^>l  V X2>1  V x3  > 1 V x4  > 1 . 

If  we  apply  Theorem  3.1  without  strengthening  and  choose  0*  = (1  - a^) 
i = 1,2, 3, 4,  we  obtain  the  cut 

2 ,2  1 ^ . 

3 X5  5 X6  3 x7  - 1 * 


whose  j-th  coefficient  is 


-a 


Of  j 


max 


1L 


ie{l, 2,3,4}  1 ” ai0 


To  apply  the  strengthening  procedure,  we  note  that  each  Xj,  j ■ 1,2, 3, 4 
is  bounded  below  by  0.  Using  a*  = 1 (which  satisfies  (6.11)  since 


o1(3q  - b*)  = oi[l  - aiQ  - (-a1Q)]  « 1,  1 ■ 1,2, 3, 4,  and  0*3* 


° (1  'ai0)  > °* 


we  obtain 
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min  max 
meM 


min  max 
meM 


t ! < I + ”!>■  ! <’  l + "21-  ! <_  I + ”3>f  5 ('  f + ”«)! 


t I (-I  + V-  ! < 6 + V-  3 ( i + ”3)>  I ( 6 +”4,) 


6 , 1 


meM 


Next  we  apply  the  above  Algorithm  for  calculating 


For  j = 5:  y = - ^ ml  = 34  * m2  “ " 102  ’ m3  51  ’ “A  51 


* 32 

m„  = 


* 11 
m.  = 


Thus  our  starting  values  are  [m^]  - 0,  [n^]  - ”*»  Im3^  ^A^ 

Since  k = -(-1)  - (~1)  = 2>  the  Iterative  step  is  applied  twice: 


i min  I - - - - 1 = - 7 , s = 2;  m = 0,  m2  = -1+  1 - 

1.  min  l 4»  3»  5 -I  4 ’ 1 *■ 

m,  = -1,  m.  = 0. 


2.  min 


r I 5.  2 2 1 = - I s - 1;  m.  • 1,  m,  • 0,  m = -1,  ® - °* 

1“  5>  4»  3*  5 ■>  5 1 1 J 


These  are  the  optimal  m^  and 

r 1 I _ A . A = _ I 

P5  - max  [-  5,  - 4,  3»  5 J 5 

For  J - 6;  Y - - [<1  - -1.  ^ ^ = °’  ^A1  = «;  k 2. 

1.  min  ( f,  » J.  f.  5 ) " " 1 S = 25  "l  " ’li  ”2  = ra3  = °'  “A  = ° 

, , f 2 521,1  . 4;  m = -1,  m - o,  m = 0,  m -1. 

2.  min  l J>  4»  3'  5 J 5 * 

, A i .A  1,1. 

P6  = max  5.  - 4>  3»  5 3 5 


2 ^ 

— For  J - 7:  Y * “ yj>  ^1^ 


= 0,  tm*]  = -1,  K1  = “l»  K1  = ”l;  k = 3‘ 


, ,„r  I I I I 1 = i s = 1;  m.  * 1,  m, 

1.  min  l5»4»3»5J  5 1 L 


■1,  m_  = -1»  m.  «*  ”li 
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„ , r 7 1 1 1 , 1 , 

2.  min  l 4'  3*  ^ 1 c 5*  6 a ‘‘i  »j 


, f 7 1 1 7 , 1 

3.  min  [ ^ "7  > 5 J "4*  3 ~ ulj 

, 1 1 5 1,1 

P;  - max  t 5'  41  ‘ 3’  5 ' “ 4 ' 


L,  m2  = -1,  “ -1,  » 0; 


1,  m2  “ n3  “ “1»  m4  " 0* 


Thus  the  strengthened  cut  is 


1 ^ 1 . 1 ^ , 

T Xc  + T X,  + 7 X-,  > 1 

5 5 5 6 4 7 — 


The  frequently  occurring  situation , when  jQj  - 2,  deserves  special 
mention.  In  this  case  the  coefficients  |3  , .1  e Jj,  are  given  t 


(7.15) 
where 

(7.16) 


3rC  7, 

= min  + <mf)>),  \,^c'2  " ^ » 


m. 


2 * 1°  I 
xi>>;2 


with  , Xt>  i = 1,2,  defined  by  (7.12),  and  <ji1q>  “ the  Sinallest  integer  >m 
The  optimal  value  of  = -m^  is  either  <m.Q>  or  tn^] , according  to  whether 
the  minimum  in  (7.15)  is  attained  for  the  first  or  the  second  term. 

The  strengthening  procedure  discussed  in  this  section  produces  the 
seemingly  paradoxical  situation  that  weakening  a disjunction  by  adding  a new 
term  to  it,  may  result  in  a strengthening  of  the  cut  derived  from  the 
disjunction;  or,  conversely,  dropping  a term  from  a disjunction  may  lead  to 
a weakening  of  the  inequality  derived  from  the  disjunction.  For  instance, 
if  the  disjunction  used  in  Example  7.2  is  replaced  by  the  stronger  one 

Xj  > 1 V x2  > 1 V > 1 , 

then  the  cut  obtained  by  the  strengthening  procedure  is 

1 x 2 x 1 ^ , 

“ 5 x5  + 5 x6  + 4 x7  ^ 1 • 


,Q. 
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which  is  weaker  than  the  cut  of  the  example,  since  the  coefficient  of  x, 

6 

is  2/5  instead  of  1/5.  The  explanation  of  this  strange  phenomenon  is  to 
be  sought  in  the  fact  that  the  strengthening  procedure  uses  the  lower  bounds 
on  each  term  of  the  disjunction.  In  Example  7.2,  besides  the  disjunction 

> 1 V X2>1  V x^  > 1 V > 1,  the  procedure  also  uses  the  information 
that  x^  > 0,  i « 1,2, 3, 4.  When  the  above  disjunction  is  strengthened  by 
omitting  the  term  x^  > 1,  then  the  procedure  does  not  any  more  use  the  infor- 
mation that  x,  >0. 

4 — 

8.  Some  Frequently  Occurring  Disjunctions 

As  mentioned  earlier,  one  of  the  main  advantages  of  the  disjunctive 
programming  approach  is  that  it  can  make  full  use  of  the  special  structure 
inherent  in  many  combinatorial  problems.  In  [ 8 ],  [ 9 ] cutting  planes  are 
derived  from  the  logical  conditions  of  the  set  partitioning  problem,  the 
linear  complementarity  problem,  the  two  forms  of  representation  for  nonconvex 
separable  programs,  etc.  More  general  complementarity  problems  are  discussed 
in  [33],  [34],  Here  we  illustrate  the  procedure  on  the  frequently  occurring 
logical  condition  (where  x^  > 0 integer,  i e Q) 

(8.1)  £ x - i, 

ieQ 

often  called  a multiple  choice  constraint. 

If  all  the  problem  constraints  are  of  this  form,  we  have  a set  par- 
titioning problem.  But  the  cut  that  we  derive  uses  on’y  one  condition 

(8.1) ,  so  it  applies  to  arbitrary  integer  programs  with  at  least  one 
equation  (8.1).  It  also  applies  to  mixed-integer  programs,  provided  the 
variables  x^,  ieQ,  are  integer-constrained. 


Let  I and  J be  the  index  sets  for  the  basic  and  nonbasic  variables  in 


a basic  feasible  noninteger  solution  of  the  form 


(8.2) 


ai0  + E aij(“Xj)  * 1 * 1 • 

J ^ J 


In  [8],  [9],  several  cutting  planes  are  derived  from  the  disjunction 


Xi  = 1 


ieQ  \ x^  = °>  Vh  e Q \ [i  } 


Here  we  describe  another  cut,  which  in  most  cases  turns  out  to  be 
stronger  than  those  mentioned  above.  It  is  derived  from  the  disjunction 


(8.3) 


£xt=OV  E = 0 , 


clearly  valid  for  any  partition  Q sense  of  being  satisfied 

by  every  integer  x satisfying  (8.11. 


Denoting 


(8.4) 


ielHQ 


E . t - 1,2;  , 

rrw  J 


(8.3)  can  be  written  as 


£ 6^*3  > ej 

jeJ  J J 


E x = 0 
jcJ(1Q1  J 


jej  ^ 


E x.  = 0 


which  implies  the  disjunction 


(8.5) 


J 

jeAc^  J J 


2 „2 

V E p.x,  > , 

jeAQ2  3 3 


1 
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ir 

with  pQ  > 0,  k = 1,2.  Note  that  once  the  sets  I 0Qk>  k = 1,2,  are  chosen 

the  sets  J flQ^  and  J n can  be  "optimized,"  in  the  sense  of  putting  an 

index  j eJOQ  into  if  (p^/p^)  > (p^/p^),  and  in  Q2  otherwise.  Using 

k k 

this  device  while  applying  Theorem  3.1  to  (8.5)  with  multipliers  6 ■=  1/p^  , 

k - 1,2,  we  obtain  the  cut 

(8.6)  E PjXj  > 1 , 

with  coefficients 


We  now  apply  the  strengthening  procedure  of  section  7 to  the  coefficients 
Pj»  J\Q  (the  coefficients  indexed  by  J f|Q  can  usually  not  be  further 
strengthened).  This  of  course  assumes  that  all  x.,  j eJ\Q,  are  integer 
constrained.  A lower  bound  on  E p.x.  is  p - 1,  for  k = 1,2,  since 

jeAQk  J J ° • 


' Po  + E pk(-x.)  < 1 • k “ 1»2- 

ieinQk  jeJ\Qk  J J 


E x. 


The  multipliers  o = 1,  k = 1,2,  satisfy  condition  (7.11)  of  Corollary  7.1.1, 


since 


okIPo  - (pk  - 1)]  = 1,  akpk  >0,  k = 1,2, 


and  thus  the  j-th  coefficient  of  the  strengthened  cut  becomes  (Corollary  7.1.1) 
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\ 


P 


i 


min  max 
meM  ke{l,2} 


Jii 

N k 


l 


P, 


0 


with  M defined  by  (7.5),  with  | Q | =2.  Applying  the  closed  form  solution  (7.15) 
to  the  minimization  problem  involved  in  calculating  p (in  the  special  case  of  a 
disjunction  with  only  two  terms),  we  obtain 


kDk  ok  -v  / k k. -1  1 , , o 

c*l  * O Pj  = Pj  , “ (<*  Pq)  = "k  ’ k c;]L»2» 

and  hence 


where 

(8.8) 


min 


+ <m 


* 

0 


> 


1 

0 


* 

m0 


fiV  q1q2 

Yo  " P^Pp 
1 , 2 
Pq  + Pq 


We  have  thus  proved 

Theorem  8.1.  If  (8.2)'  is  a basic  feasible  noninteger  solution  of  the 


linear  program  associated  with  an  integer  program  whose  variables  have  to 
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% 


where  the  k = 1,2,  j e J,  are  defined  by  (8.4)  and  m^ 
We  illustrate  this  cut  on  a set  partitioning  problem 
special  case  of  the  Theorem. 


is  given  by  (8.8) . 
, which  is  a 


Examp le  8.1.  Consider  the  set  partitioning  problem  whose  cost  vector 
is  c = (5 ,4,3 ,2 ,2 ,3 , 1, 1, 1,0)  and  whose  coefficient  matrix  is  given  in  Table  8.1, 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1 

1 

1 

1 

1 

2 

1 

1 

1 

3 

1 

1 

1 

1 

1 

4 

1 

1 

1 

1 

5 

1 

1 

1 

1 

1 

1 

Table  8.1 


1 2 

The  linear  programming  optimum  is  obtained  for  x.  = xc  = xn  = -r,  x,  = — . 

4 o j J 7 3 

xio  = l>  and  Xj  = 0 for  all  other  j.  The  associated  system  (8.2)  is  shown 
in  the  form  of  a simplex  tableau  in  Table  8.2  (artificial  variables  have  been 
removed ) . 


x 


0 


x 


10 


x 


8 


x 


7 


x 


9 


1 

'X1 

“x6 

“x3 

'X5 

,-x2 

-2 

5 

2 

1 

1 

3 

1 

1 

0 

0 

0 

1 

1 

1 

2 

2 

1 

1 

3 

3 

3 

3 

3 

3 

2 

1 

1 

1 

1 

2 

3 

3 

3 

3 

"3 

3 

1 

1 

2 

2 

1 

1 

3 

3 

3 

3 

' 3 

3 

1 

2 

2 

1 

1 

2 

3 

3 

3 

" 3 

3 

3 
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We  choose  the  disjunction  corresponding  to  row  5 of  the  matrix  shown  in 
Table  8.1,  which  is  of  the  form  (8.1),  with  Q = {3 ,4 , 5 , 6 , 8 , 9 } . 

We  define  I 0Q1  - {4,8},  I PQ2  = {9},  and  we  have 


J = 1 j =6  j = 3 j ” 5 j=2  j -0 


31 

J 3 


R2  1 

pj  3 


Since  J 0Q  = {3,5,6},  we  need  the  values  for  jeJ\Q  = {1,2}. 

* 2 * 2 
They  are  m0(l)  = ^ > mQ(2)  = ~ 3 * Hence 


2 2 1 
•|+1  3-°l  _ 1 


= min  \ 0 , max 


T-Tf 

l 3 3 J 


min  \ 0 , max 


which  is  considerably  stronger  than  the  traditional  cuts  that  one  can 
derive  from  Table  8.2,  and  it  actually  implies  that  the  nonbasic  variable  x2 

has  to  be  1 in  any  integer  solution. 

Dual  cutting  plane  methods  have  been  found  reasonably  successful  on 
set  partitioning  problems.  Using  stronger  cuts  can  only  enhance  the 
efficiency  of  such  methods,  since  the  computational  cost  of  the  cut  (8.9) 
is  quite  modest. 

• * . ' " •» 
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9.  Combining  Cutting  Planes  with  Branch  and  Bound 

The  disjunctive  programming  approach  offers  various  ways  of  combining 
branch  and  bound  with  cutting  planes,  some  of  which  are  currently  the 
object  of  computational  testing.  Here  we  discuss  one  feature  which  seems 
to  us  crucial. 

For  any  problem  P,  let  v(P)  denote  the  value  of  P (i.e.,  of  an  optimal 
solution  to  P). 

Suppose  we  are  using  branch  and  bound  to  solve  a mixed-integer  0-1 

program  P,  stated  as  a maximization  problem.  If  { P^} ^ gp  is  the  set  of 

active  subproblems  (active  nodes  of  the  search  tree)  at  a given  stage  of 

the  procedure  and  v(P^)  is  the  available  upper  bound  on  v(P^)  (assume,  for 

the  sake  of  simplicity,  that  v(P^)  = v(U\),  where  1^  is  the  linear  program 

associated  with  P.),  then  max  v(?i>  is  an  upper  bound  on  v(P) . Also, 

1 ieQ 

v(P) , the  value  of  the  best  integer  solution  found  up  to  the  stage  we  are 
considering,  is  of  course  a lower  bound  on  v(P);  i.e.,  at  any  stage  in  the 
procedure  , 

(9.1)  v(P)  < v(P)  < max  v(P. ) . 

ieQ  1 

Hence  the  importance  of  finding  good  bounds  for  both  sides  of  (9.1), 

It  is  a crucial  feature  of  the  approach  reviewed  here  that  it  can  be 

used  to  derive  a cutting  plane  from  the  optimal  simplex  tableaus  associated 

with  the  subproblems  LP^ , ieQ,  which  provides  an  upper  bound  on  v(P)  at  least 

as  good  as,  and  often  better  than  max  v(P.). 

ieQ  1 

Let  the  linear  program  LP  associated  with  the  mixed-integer  0-1  program 
P have  an  optimal  solution  of  the  form 

*h  = aho + ^ 
jeJ 


(9.2) 


VV 


h e 1 U {0} 
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where  I and  J are  the  index  sets  for  the  basic  and  nonbasic  variables 
respectively,  and  let  I,  and  be  the  respective  index  sets  for  the  integer 
constrained  variables.  Here  a ^ > 0,  he  I,  and  < 1,  h e 1^ . Further, 

since  P is  a maximization  problem  and  the  solution  (9.2)  is  optimal,  a > 0, 

J e J. 


Now  let  (PjJkeQ  ^e  C^e  set  active  subproblems,  and  for  keQ,  let  the  optimal 
solution  to  LP^,  the  linear  program  associated  with  P^,  be  of  the  form 


(9-3), 


*h  = ah0+  .^k  ahj("xj)  * 


h cl  U{0}  , 


If  V 

where  I , J are  defined  with  respect  to  LP^  the  same  way  as  I,  J 

with  respect  to  LP.  Again  aQ  > 0,  Vj  eJ  , since  each  LPk  is  a maximization 

problem. 

In  order  to  derive  a valid  cutting  plane  from  (9.3)^,  keQ,  we  view 
the  branching  process  as  the  imposition  of  the  disjunction 


(9.4) 


V 

keQ 


where  Ax  > b stands  for  the  system 


ahj*j  2 %o  - 1 • 


h e I 

h € 1^  y 


expressing  the  conditions  > 0,  h € I,  x^  < 1,  he  1^,  while  each 
k k 

D x > dg  is  composed  o i inequalities  of  the  form 


with  Ax  > b,  x > 0,  defines  P^. 

Now  consider  the  cut  derived  from  (9.4)  on  the  basis  of  Theorem  3.1, 

with  the  optimal  dual  variables  (obtained  by  solving  LP^)  used  as  the 

lc  lc  1c 

multipliers  9 , k cQ.  If  for  kcQ  we  dencte  by  (u  ,v  ) the  optimal  dual 

vector  associated  with  the  k-th  term  of  (9.4  "»id 
(9.5)  ak  = ukDk  + vkA  , 


k kjk  ■ Si  , 


a0  = u d„  + v 


0 


and  if  > 0,  k e Q , then  according  to  Theorem  3.1,  the  inequality 

k 


(9.6) 


E (max  — f ) x . > 1 

, _ , _ k i — 

jcJ  keQ  qtq 


is  satisfied  by  every  x satisfying  (9.4),  i.e.,  by  every  feasible  integer 

solution.  The  condition  0q  > 0,  keQ,  amounts  to  requiring  that  v(LP^)  < v(LP), 

k k 

i.e.,  that  the  "branching  constraints"  D x > d^  force  v(LP  ) strictly  below 

v(LP),  VkeQ.  This  is  a necessary  and  sufficient  condition  for  the  procedure 

discussed  here  to  be  applicable.  Should  the  condition  not  be  satisfied  for 

some  keQ,  one  can  use  a different  objective  function  for  LP^  than  for  the 

rest  of  the  subproblems  — but  we  omit  discussing  this  case  here.  Note  that, 

since  b < 0,  v^  <0,  Vk  eQ,  and  thus  Qk  > 0 implies  ukdk  >0,  Vk  *Q. 

k k 

Since  the  multipliers  (u  ,v  ) are  optimal  solutions  to  the  linear 
programs  dual  to  LP^,  keQ,  they  maximize  the  right  hand  side  coefficient 
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k k k 

Oq  of  each  inequality  a x > q-q  underlying  the  cut  (9.6)  subject  to  the 
condition  that  < a0j,  VjeJ. 

We  now  proceed  to  strengthen  the  inequality  (9.6)  via  the  procedure  of 
section  7.  To  do  this,  we  have  to  derive  lower  bounds  on  a x,  k eQ.  We 
have 

k k k k ,k.  , k . . . 

a x - a0  = u (D  x - dQ)  + v (Ax  - b) 

„ k _ k ,k.  . k 

> u (D  x - dg)  > -u  e 

where  e = (1,...,1).  The  first  inequality  holds  since  Ax  - b > 0 for  all  x 

satisfying  (9.4),  while  the  second  one  follows  from  the  fact  that  each 

Ic  k 

inequality  of  the  system  D x • d^  > 0 is  either  of  the  form  -x^  > 0 or  of  the 
form  - 1 > 0,  and  in  both  cases  -1  is  a lower  bound  on  the  value  of  the 
left  hand  side.  Thus 

(9.7)  a^x  > q-q  - u^e  , k eQ 

k k k 

holds  for  every  x satisfying  (9.4).  Note  that  u d^  > 0 implies  u ^0  and 
k k 

hence  u > 0 implies  u e > 0,  k eQ. 

We  now  apply  Corollary  7 . 1 . 1 to  the  system 

k k k , 

x > - u e , keQ  , 

(9.8)  x > 0 , 

v (a*  > c£>  , 

kcQ  U 

Xj  integer,  jeJ^. 

k k 

We  choose  a = 1/u  e,  keQ,  which  satisfies  condition  (7.11)  of  the 
Corollary: 
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k r k , k k 


(l/u^e)  [qtq  - (o“  - u"e)]  = 1 , (l/uke)akfl  > 0 


The  strengthened  cut  is  then 


(9.9) 


(9.10) 


£ P«X4  > 1 * 

jcJ  J J 


(a^/u^e)  + m. 


min  max  ■— - — - — 

mcM  keQ  c*q/u  e 


k » 
keQ  a0 


j e J\J. 


where 


(9.11) 


M = {m  c R ' | E m^  > 0,  m^  integer,  keQ]  . 
keQ 


The  values  of  a j > Oq  and  u e needed  for  computing  the  cut  coefficients,  are 
readily  available  from  the  cost  row  of  the  simplex  tableaux  associated  with 
the  optimal  solutions  to  LP  and  LP^,  keQ.  If  the  latter  are  represented  in  the 
form  (9.2)  and  (9.3)^  respectively,  and  if  d^  and  a^  denote  the  j-th  column 

It 

of  D and  A of  (9.4),  while  S,  is  the  row  index  set  for  D , keQ,  we  have 


for  all  keQ 


„ k k k 
a0j  - u dj  - v a . 


aoj  " aj  ’ 


0+"J  . 


j c jk  n j 


J ej  nsk  , 


since  the  indices  j e correspond  to  the  slack  variables  of  the  system 
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k k 

-D  x < -d  , whose  costs  are  0 (note  that  S,  OJ  = 0 by  definition).  Further, 

“ U K 

k k 

for  j e J\J  =JflI  the  reduced  cost  of  x^  in  LP^  is  0,  hence  for  all  k cQ 

n k,k  k 

0 ■ .0J  - » dj  - v as 


Finally, 


00 


a0j  ' aj  • 


k k k, 
a„„  — u d - v d 
00  0 


00 


k 

*0  ’ 


Vk  e Q . 


From  the  above  expressions  we  then  have  for  keQ, 


(9.12) 


Of, 


“oj  ‘ aoj  ' 3ejnjk 


Oj  ’ 


j e J\J 


k 

er„  = a 


00  00  ’ 


and 


(9.13) 


k 

u e = 


E a 

ieJkns, 


Oi 


(since  uk  = 0,  VieS^fll^). 

Hie  representation  (9.3)^  of  the  optimal  solution  to  LP^  assumes  that 
the  slack  variable  j e of  each  "branching  constraint"  x^  < 0 or  x^  > 1 
that  is  tight  at  the  optimum,  is  among  the  nonbasic  variables  with  a^  > 0. 
If  one  prefers  instead  to  replace  these  slacks  with  the  corresponding 


structural  variables  x^  and  regard  the  latter  as  "fixed"  at  0 or  1,  and  if 
F^  denotes  the  index  set  of  the  variables  fixed  in  LP^,  the  reduced  costs 
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k k 

a^,  i eJ  flF^  are  then  the  same,  except  for  their  signs,  as  agj>  j cJ  flS^, 

and  the  only  change  required  in  the  expressions  derived  above  is  to  replace 

(9.13)  by 


(9. 13') 


l |a0i' 

ieJ  TFk 


Of  course,  in  order  to  calculate  a cut  of  the  type  discussed  here,  one 
needs  the  reduced  costs  a^  for  both  the  free  and  the  fixed  variables. 

We  have  thus  proved  the  following  result. 

Theorem  9.1.  If  LP  and  LP^,  k cQ,  have  optimal  solutions  of  the  form  (9.2) 
and  (9.3)^  respectively,  with  a ^ > a^,  keQ,  then  every  feasible  integer 
solution  satisfies  the  inequality  (9.9),  with  coefficients  defined  by  (9.10), 
(9,11),  (9.12)  and  (9.13)  [or  (9.13')]. 

In  the  special  case  when  |q|  = 2 and  LP^,  LP^  are  obtained  from  LP  by 
imposing  x^  < 0 and  x^  > 1 respectively  (for  some  i £ I such  that  0 <  2  3i0  < 1 )» 
the  definition  of  (3^  for  j e becomes 


(9 . 10 /) 


f (<*W) 


Pj  = min 


i1>  + <^>  c^-^1 
1/1  » 2,  2 
Qf0/ui  or0/u.  j 


(9.14) 


2 1 12 
*1*0  ' Q1a0 


2 1 12  * 
uia0  + uia0 


We  now  state  the  property  of  the  cut  (9.9)  mentioned  at  the  beginning 
of  this  section. 

Corollary  9.1.1.  Adding  the  cut  (9.9)  to  the  constraints  of  LP  and 


performing  one  dual  simplex  pivot  in  the  cut  row  reduces  the  value  of  v(LP) 
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from  aQQ  to  aQ0  such  that 


(9.15) 


“00  S “ “00 


Proof.  For  each  keQ,  a^  = aQQ  - Oq  • Now  suppose  (9.1)  is  false, 
*-e-’  i00>a00’  Tk,Q*  ^ 


a„„  = a. 


00  00  jej|j  >0 


”,ln  “oi/ej  > aoo  ' 'Jo  ■ vk“Q  • 

t I q J J 

j' 


and  hence  for  all  keQ, 


(9.16) 


where 


(9.17) 


“o  > m,in^  aoj/pj 

jeJ|Pj>0  J 


8 S 

> min  aQj  <a0/aj  ^ * 


jeJ |or j > 0 


s 

“o 


h 

Ofj 


= max 


heQ  a. 


0 


The  second  inequality  in  (9. 16) holds  since  the  cut  (9.9)  is  a 
strengthened  version  of  (9.6),  in  the  sense  that 


Pj  < max  , 


j e J 


heQ  q-q 


Now  suppose  the  minimum  in  the  second  inequality  of  (9.16)  is  attained 
for  j = t.  Since  (9.16)  holds  for  all  keQ,  we  then  have 

«b  > aot(Q'o/at) 
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s s s 

or  (since  >0,  <*g  > 0),  ot  > agt*  But  this  contradicts  the  relation 
a®  < aQt  implied  by  (9.12).  || 

Note  that  the  Corollary  remains  true  if  the  strengthened  cut  (9.9) 
is  replaced  by  its  weaker  counterpart  (9.6).  in  that  case,  however,  (9. 15) 
holds  with  equality.  The  remarkable  fact  about  the  property  stated  in  the 
Corollary  is  that  by  using  the  strengthened  cut  (9.9)  one  often  has  (9.15) 
satisfied  as  strict  inequality.  More  precisely,  we  have  the  following 

Remark  9.1.  (9.15)  holds  as  strict  inequality  if  < cr^/og  > where  s 

is  defined  by  (9.17),  and  t by 


a0tfa0/o£)  = 


mm 
s 


jeJ joj  > 0 


Note  that  for  (9.15)  to  hold  as  strict  inequality  the  pivot  discussed 
in  the  Corollary  need  not  occur  on  a "strengthened"  cut  coefficient.  All 
that  is  needed,  is  that  the  coefficient  on  which  the  pivot  would  occur  in 
case  the  unstrengthened  cut  were  used,  should  be  "strengthened"  (i.e., 

reduced  by  the  strengthening  procedure). 

The  significance  of  the  cut  of  Theorem  9.1  is  that  it  concentrates 
in  the  form  of  a single  inequality  much  of  the  information  generated  by 
the  branch  and  bound  procedure  up  to  the  point  where  it  is  derived,  and 
thus  makes  it  possible,  if  one  so  wishes,  to  start  a new  tree  search  while 
preserving  a good  deal  of  information  about  the  earlier  one. 
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Theorem  9.1  and  Its  Corollary  are  stated  for  (pure  or  mixed)  0-1  programs; 

but  the  0-1  property  (as  opposed  to  integrality)  is  only  used  for  the 

k k 

derivation  of  the  lower  bounds  on  D x > d^,  keQ;  hence  it  only  involves  those 
variables  on  which  branching  has  occurred.  The  results  are  therefore  valid 
for  mixed-integer  programs  with  some  0-1  variables,  provided  the  strengthening 
procedure  is  only  used  on  cuts  derived  from  branching  on  0-1  variables. 


Example  9.1.  Consider  the  problem  in  the  variables  x^  > 0,  j = 1,...,6; 
Xj  = 0 or  1,  j = 1,4;  x^  integer,  j = 2,3,  whose  linear  programming  relaxation 
has  the  optimal  solution  shown  in  Table  8.1. 


1 

-X3 

“x4 

-X5 

"X6 

1.1 

2.0 

.2 

.05 

1.17 

.2 

.4 

1.3 

-.01 

.07 

.9 

-.3 

.4 

-.04 

.1 

Table  9.1 


If  we  solve  the  problem  by  branch  and  bound,  using  the  rules  of  always 
selecting  for  branching  (a)  LP^  with  the  largest  aQ0,  and  (b)  xi  with  the 
largest  max  {up  penalty,  down  penalty],  we  generate  the  search  tree  shown 
in  Fig.  9.1.  The  optimal  solution  is  x = (0,2,1,0,20,0),  with  value  -1.9, 
found  at  node  6.  To  prove  optimality,  we  had  to  generate  two  more  nodes, 
i.e.,  the  total  number  of  nodes  generated  (apart  from  the  starting  node)  is  8. 


Suppose  now  that  after  generating  the  first  four  nodes,  we  wish  to  use 


the  available  information  to  derive  a cut.  At  that  point  there  are  three 


active  nodes,  associated  with  LP^,  for  k = 2,3,4.  The  corresponding  reduced 
k k 

cost  coefficients  a^^ , jeJ  , are  shown  in  Table  9.2.  The  slack  variables 
of  the  "branching  constraints"  x^  < 0,  x^  > 1,  x^  < 0,  and  x^  > 1 are 
denoted  by  x,,  x„. 


xQ  and  xin  respectively 


and  the  cost 


k k k 

The  coefficients  Oj.  <*o  and  u e,  extracted  from  Table  9.2 
row  of  Table  9.1,  are  as  follows: 


. - 2 2.2  2 2 
^1-  “o  = a00  - a00  » 4;  u e = = aQ8 


2 2 
5;  a3  = a03  a03  " ”2, 


a4  K a04  ” a04  ”6,5,  a5  = a05  = ,05»  “6  " a06  ” a06  = ",35: 

3333  33  3 3 

k = 3.  Qg  = l‘>  u e = u.,  + = 11.3;  a3  = 2,  = .2,  - -.05,  a ^ * .35; 

4 4 4 4 4 4 4 

K=  4:  aQ  = 5.7;  u e = a1Q  = 6.7;  a3  = -2,  = .2,  a5  “ .05,  a6  - -.35. 

The  coefficients  of  the  strengthened  cut  are  shown  below,  where 
M = {m  e R | m^  + m2  + m^  > 0,  m^  integer,  i = 1,2,3}. 


(3.  **  min  max 

J meM 


{ 


(-2/5)  + mx  (2/11.3)  + m2  (-2/6.7)  + m3 

4/5  * 1/11.3  ’ 5. 7/6. 7 / = *75» 


with  m = (1,-1,0). 


P/.  " 


/ 


min  max  -j 
meM 


4/5 

with  m = (1,-1,0) 


(-6.5/5)  + mx  (.2/11.3)  + m2  (.2/6.7)  + m , 

475  ’ 1/11.3  » 5. 7/6. 7 j = ,035 


f .05  -.05  .05  \ 

max  l 4 * 1 ’ 5.7  J 


= .0125 


Pa  = 


/ -.35  .35  -.35 

max  1 “4”  * T * T77  S 


= .35 


Adding  to  the  optimal  tableau  of  LP  the  cut 

.75  x^  + .035  x^  + .0125  x,.  + .35  x^  > 1 


produces  Table  9.3  and  the  two  pivots  shown  in  Tables  9.3  and  9.4  produce 
the  optimal  Tableau  9.5.  Thus  no  further  branching  is  required. 


-xc  -x. 


1.1 

2.0 

.2 

.05 

1.17 

.2 

.4 

1.3 

-.01 

.07 

.9 

-.3 

.4 

-.04 

.1 

-1.0 

lv.75 

-.035 

-.0125 

-.35 

Table  9.3 


-1.57 

2.67 

.106 

.01675 

1.117 

-.333 

.533 

1.281 

-.016751 

-.117 

1.3 

-.4 

.414 

-.035 

.240 

1.33 

-1.33 

.047 

.01675 

.467 

Table  9.4 


-1.9  2.687 

1.387  1.0 

1.0 

20.0  -31.8 

-73.0  -60.0 

7.0 

2.0 

-2.0 

1.0 

1.0 

Table  9.5 
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10.  Disjunctions  from  Conditional  Bounds 


In  solving  pure  or  mixed  integer  programs  by  branch  and  bound,  the 
most  widely  used  rule  for  breaking  up  the  feasible  set  Is  to  choose  an 
integer-constrained  variable  x^  whose  value  a at  the  linear  programming 
optimum  is  noninteger,  and  to  impose  the  disjunction  (x^  < [alol)V(x1  > [atQ]  +1). 
It  has  been  observed,  however,  that  in  the  presence  of  multiple  choice 


constraints,  i.e.,  of  constraints  of  the  form 


E - 1 . 


it  is  more  efficient  to  use  a disjunction  of  the  form 


( S x£  = 0)  V ( E x.  = 0)  , 


where  U Q2  = Q,  0 Q2  = 0,  and  and  Q2  are  about  equal  in  size. 

This  is  just  one  example  of  a situation  where  it  is  possible  to  branch 
so  as  to  fix  the  values  of  several  variables  on  each  branch.  The  circumstance 
that  makes  this  possible  in  the  above  instance  is  the  presence  of  the 
rather  tight  multiple  choice  constraint.  More  generally,  a tightly 
constrained  feasible  set  makes  it  possible  to  derive  disjunctions  stronger 
than  the  usual  dichotomy  on  a single  variable.  On  the  other  hand,  the  feasible 
set  of  any  integer  program  becomes  more  or  less  tightly  constrained  after 
the  discovery  of  a "good"  solution  (in  particular,  of  an  optimal  solution), 
provided  that  one  restricts  it  to  those  solutions  better  than  the  current 
best.  Such  a "tightly  constrained"  state  of  the  feasible  set  can  be 
expressed  in  the  form  of  an  inequality 


TTX  < TT 
— O 
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with  tt  > 0,  and  with  a relatively  small  ttq  > 0.  One  way  of  doing  this, 
if  the  problem  is  of  the  form 

(P)  min{cx|Ax  > b,  x > 0,  x^  integer,  jeN}  , 

with  c integer,  and  if  z^  is  the  value  of  the  current  best  integer  solution, 
ia  to  find  a set  of  multipliers  u such  that 

uA  < c , u > 0 , 

and  define 


tt  = c-  uA,tt  = z„  - ub  - 1. 

* o U 

Then  multiplying  Ax  > b by  -u  and  adding  the  resulting  inequality, 

-uAx  < - ub,  to  cx  < Zy  - 1,  yields  the  inequality 

TTX  < TT  , 

— O 

satisfied  by  every  feasible  integer  x such  that  cx  < z^.  Here  n > 0,  Tio  > 0, 
and  the  size  of  ttq  depends  on  the  gap  between  the  upper  bound  z^  and  the 
lower  bound  ub  on  the  value  of  (P) . 

Now  suppose  we  have  such  an  inequality  nx  < tt  . Without  loss  of 
generality,  we  may  assume  that  > 0,  V-  j (by  simply  deleting  the  zero 
components).  Then  the  following  statement  holds.  [12], 

Theorem  10.1.  Let  tt  e Rn,  ttqcR,  (tt,  tt  ) > 0,  N = {l n},  and  for 

i = l,...,p,  1 < P < n,  let  Q^CN,  Q t 0,  with 


"jCD 


= min  n 

jeQi  J 


If  the  sets  Q^,  i = 1 


(10.1) 


Z TT 


j(i)  - n 


j 


satisfy  the  conditions 
, jeN, 


and 


(10.2) 


P 


E TT 
i=l 


Jd) 


> TT 

O 


» 


then  every  Integer  vector  x > 0 such  that  nx  < tt^,  satisfies  the  disjunction 


(10.3) 


V (x  = 0,  JeQ  ) 
1=1  J 


Proof.  Every  integer  x > 0 which  violates  (10.3)  satisfies 


(10.4) 


E x,  > 1 > 1 - 1 , . . . ,p. 

jeQi  3 


4.L 

Multiplying  by  the  1 Inequality  of  (10.4)  and  adding  up  the 

resulting  inequalities,  one  obtains 

P 


(10.5) 


Further, 


Jc(^Q  i'JeQi 
i=l 


(2_I  "j(1)>*j  > ^Jd)  • 


E n x > n.x. 

JeN  J J Jc0Q  J J 
i=l 


[£t°"  <10•I,, 


P 

> En,.,.  [from  (10.5)] 

i=l 


> n [from  (10.2) ] . 

o 


Q.E.D. 

One  way  of  looking  at  Theorem  10.1  is  as  follows.  Suppose  the  constraints 
of  an  integer  program  which  include  the  inequality  irx  < ttq,  were  amended 
by  the  additional  constraints  (10.4). 
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From  the  proof  of  the  Theorem,  these  inequalities  give  a lower 
bound  on  ttx  which  exceeds  ttq;  this  contradiction  then  produces  the  dis- 
junction (10.3).  Since  the  inequalities  (10.4)  are  not  actually  part  of  the 
problem,  we  call  the  bound  on  nx  derived  from  them  a conditional  bound, 
and  the  disjunction  obtained  from  such  bounds,  a disjunction  from  conditional 
bounds . 

Example  10.1.  The  inequality 

9x1+8x2+8x3+7x4+7x5+6x6+6x7+5xg+5x9+5x10+4x11+4x12+3x13+3x14H-3x15+2x16 

+2x17  < 10, 

together  with  the  condition  x > 0,  x^  integer  ¥ j,  implies  the  disjunction 
(Xj  = 0,  j = 1,2, 3,4, 5, 6, 7)  V (Xj  - 0,  j * 1,8,9,10,11,12,13,14)  V 

v (Xj  = 0,  j « 2,3,8,9,10,15,16,17). 

Indeed , 

"jU)  = mlnt9»8*8*7>7»6»63  * 6» 

nj(2)  * mlnl9»5,5,5,4,4,3,3}  = 3, 

TTj  ^3^  = min{8 ,8 , 5 ,5 ,5 , 3 ,2 ,2)  = 2, 

and  (10.1),  (10.2)  are  satisfied,  since  6 + 3 + 2 > 10,  while  6 + 3 < 9 (j  « 1), 
6 + 2 < 8 (j  = 2,3),  6 < 7 ( j = 4,5),  6 < 6 (J  = 6,7),  3 + 2 < 5 ( j = 8,9,10), 
3 < 4 (j  = 11,12),  3 < 3 (j  = 13,14),  2 < 3 (j  = 15),  2 < 2 (j  = 16,17). 

Next  we  outline  a procedure  [12]  based  on  Theorem  10.1  for  systematically 
generating  disjunctions  of  the  type  (10.3)  from  an  Inequality  ttx  < ttq,  with 
TTj  > 0,  ¥ j. 


73< 


1.  Choose  some  SCN  such  that 


E tt 
JcS 


J 


> TT 

O 


but 


E TT  < TT 
j*T  3 " 


o 


for  all  TCS,  T ^ S.  Order  S * {j  (1)  , . . . , j (p)  } according  to  decreasing 
values  of  and  go  to  2. 

2.  Set 

Qx  = {jeN^  > TTj(1)] 
and  define  recursively 

Qi  " ^‘N^J  -nj(i)  + Jj,nJ(k)6JJ»  1 = 2 

k k 

where  6j  = 1 if  je^,  6j  = 0 otherwise.  The  sets  Q^,  i « 1 p, 

obtained  in  this  way,  satisfy  (9.1),  (9.2). 

In  the  above  example,  S **  {7,14,17}.  If,  on  the  other  hand,  one 
uses  S * {5,12}  (which  is  also  admissible,  since  + tt^  = 7 + 4 > 10) , 
one  obtains  the  disjunction 

(Xj  - 0.  J - 1.....5)V(X  ■=  0,  j -=  6 12). 

A disjunction  of  the  form  (10.3)  can  be  used  to  partition  the  feasible 

set  into  p subproblems,  the  kth  one  of  which  is  constrained  by 

£ l»i*lt...,k-l;  x.  * 0,  ¥ J*Q.  • 

J*Qi  3 3 k 


/ 

I 
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Another  way  of  ualng  a disjunction  of  the  form  (10.3)  is  to  derive 
cutting  planes.  This  has  been  explored  In  the  context  of  set  covering 
problems  [12]  and  found  to  yield  good  results.  In  particular,  let 
A = be  a 0-1  matrix  and  consider  the  set  covering  problem 

(SC)  min{cx|Ax  > e,  ■ 0 or  1,  JcN} 


where  e is  the  vector  of  l's  of  appropriate  dimension.  For  every  row  i 
of  A,  let 


N£  = tj«N|atj  = 1}. 

Now  suppose  a prime  cover  (a  basic  feasible  integer  solution)  x is 
known;  then  = cx  is  an  upper  bound  on  the  value  of  the  optimum.  If 
u is  any  feasible  solution  to  the  dual  of  the  linear  programming  relaxa- 
tion of  (SC),  i.e.,  any  vector  satisfying  uA  < c,  u > 0,  then  setting 
n ■»  c - uA  and  n0  = zy  " ue  one  obtains  an  inequality  ttx  < ttq  which  is 
satisfied  by  every  integer  solution  better  than  x,  and  which  can  therefore 
be  used  to  derive  a disjunction  of  form  (10.3). 

Suppose  this  is  done,  and  a disjunction  (10.3)  is  at  hand,  known 
to  be  satisfied  by  every  feasible  integer  x better  than  the  current  best. 
Then  for  each  ie[l, . . . ,'p] , one  chooses  a row  h(i)  of  A,  such  that 
is  "large"  - or,  conversely,  ^(i)^!  is  Clearly  (and  this  is 

true  for  any  choice  of  the  indices  h(i)),  the  disjunction  (10.3)  implies 


P 

(10.6)  V x > 1) 

i-l  J«»h(1,\Ql 


which  in  turn  implies  the  inequality 


Z x > 1 
JcW  J 


i 


(10.7) 


-75- 


where 


H ' .U.'WV 


i-1 


The  class  of  cutting  planes  (10.7)  obtained  in  this  way  was  shown 
in  [12]  to  include  as  a proper  subclass  the  Bellmore-Ratlif f inequalities  [16] 
derived  from  involutory  bases.  An  all  integer  algorithm  which  uses  these 
cutting  planes  was  implemented  and  tested  with  good  results  on  set  covering 
problems  with  up  to  900  variables  (see  [12]  for  details). 
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• !"T.'iT  reviews  some  recent  developments  in  the  convex  analysis  approach  ( 
Oiiri'.cr  programming.  These  developments  are  based  on  viewing  integer 

no  disjunctive  programs,  i.e.,  linear  programs  with  disjunctive 
* i ?i!tu  n,  an  approach  which  seems  to  be  fruitful  for  0-1  programming  both 
r Irn)  1 v and  practically.  On  the  theoretical  side,  it  provides  j 

tnrni  characterizations  which  offer  new  insights.  On  the  practical  j 
'•  5 ! produces  a variety  of  cutting  planes  with  desirable  properties 

‘ ' ''!  n new  ways  of  combining  cutting  planes  with  enumerative  techniques,  t 


</.)  i jinon  or  i hsv  «*  ii  onsoLrre 

'.  0 I 4 • Mill 


UNCLASSIFIED 

**c7jmiTvC^L'*7.IlVVCAT'lo51'oV~TsTs  **AOt  fWnwi  OHm  r.r.ittd) 


